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Abstract 


In  current  survivability  studies,  the  nuclear 
bursts  are  treated  as  independent  events.  Using  this 
method,  the  effect  of  thermal  radiation  from  one  fireball 
at  a  time  is  studied.  This  treatment  does  not  consider 
the  cumulative  effects  of  receiving  thermal  radiation  from 

more  than  one  fireball  at  a  time.  The  purpose  of  this 

C^’c’U  j  ~j"' 

thesis>>44h  to  develop  a  computer  program  to  model  the 

C  huc&ov  t'-cWil 

cumulative  effects' of^thermal  radiation, and  compare  these 
results  to  those  from  the  noncumulative  case. 

!fhe  scenario  studied  was  the  Peacekeeper  Dense 
Pack  missile  system.  The  missile  field  was  subjected  to  a 
walk  attack  of  2  MT  weapons  every  two  seconds,  (jpb^siming 
error  of  the  incoming  RV  was  modeled  using  a  10-cell  6 P 
cular  error  probable  (CEP)  area  around  the  designated 
ground  zero,  and  the  probability  of  damage  due  to  an  RV 
was  calculated  using  a  cumulative  log-normal  distribution 


function.  - 

'V,  ivi  1 3  3 lJj?  3  Hi rt  .  ) 

In  order  *£o  model  the  ^temperature  rise^  of  the 

missile  skin,  an  energy  balance  was  made  over  a  unit  area 


of  skin  surface  and  then  solved  using  the  thin  skin  approxi¬ 
mation  and  finite  differences .  The  resulting  equation 
gave  an  expression  for  the  skin  jtemperature  at  a  time  t 


Wv-V  -w. 


rr  .v\. ' 


(t  1‘  1 

M 

after  the  first  burst  detonated.  The  maximum  temperature 

reached  was  (thet?used  to  calculate  the  probability  of 

damage  to  the  missile  skin.  - 

To  model  cumulative  thermal  effects'^the  amount  of 

thermal  radiation  emitted  from  each  burst  was  added 

together  to  calculate  (tbe'skin  temperature.  This  method 

resulted  in  temperatures  £hat  woae^significantly  higher 

thanHthe  temperatures  calculated  for  each  burst  indepen- 

^  (iV  M  *» 

dently .  ^Bonsequently-^^umulative  thermal  effect  approved 

to  have  a  greater  region  of  no  survival  than  noncumulative 

thermal  effects  and  also  blast  effects. 
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CUMULATIVE  THERMAL  EFFECTS  IN  A 
MULTIBURST  SCENARIO 

I.  Introduction 

Background 

A  nuclear  weapon  releases  a  large  portion  of  its 
energy  in  the  form  of  electromagnetic  radiation.  This  radi¬ 
ation,  which  is  emitted  within  a  microsecond  after  the  blast, 
is  called  primary  radiation.  Since  the  primary  radiation 
is  emitted  at  tens  of  millions  of  degrees,  it  is  in  the 
soft  X-ray  region  of  the  electromagnetic  radiation  spectrum. 
Therefore,  for  a  low  altitude  or  surface  burst,  the  radia¬ 
tion  is  almost  .Immediately  absorbed  by  the  atmosphere.  As 
a  result,  the  air  is  heated  and  forms  a  fireball  that  in 
turn  emits  thermal  radiation  in  the  ultraviolet  to  infrared 
regions  of  the  electromagnetic  radiation  spectrum.  This 
secondary  thermal  radiation  travels  a  long  distance  from 
the  burst  in  a  short  time.  The  fraction  of  the  bomb  yield 
that  is  emitted  as  effective  thermal  radiation  (primary 
plus  secondary)  depends  on  the  height  of  the  burst,  the 
total  yield,  and  other  weapon  characteristics. 

As  the  thermal  radiation  travels  through  the  air, 
it  diverges  and  is  attenuated  by  absorption  and  scattering 
processes  with  air  molecules  and  other  particles.  These 
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processes  depend  on  the  wavelength  of  the  radiation  and 
atmospheric  conditions  that  vary  with  altitude  and  the  size 
and  density  of  the  interacting  particles.  Thus,  the  amount 
of  thermal  radiation  incident  on  a  target  is  determined  by 
the  total  energy  yield,  the  height  of  the  burst,  the  dis¬ 
tance  from  the  target  to  the  burst,  and  the  changing  charac¬ 
teristics  of  the  atmosphere. 

At  the  target,  a  fraction  of  the  incident  thermal 
radiation  will  be  absorbod.  For  a  given  surface  material, 
only  a  small  amount  of  the  absorbed  energy  will  be  dissi¬ 
pated  away  from  the  surface  by  conduction,  convection,  or 
re-radiation.  Therefore,  the  absorbed  energy  is  contained 
in  a  shallow  depth  of  the  target  skin,  resulting  in  high 
temperatures  that  could  damage  the  skin  material.  The  pur¬ 
pose  of  a  survivability  study  is  to  determine  what  this 
skin  temperature  will  be  in  order  to  predict  the  target's 
probability  of  survival.  If  the  skin  temperature  is  above 
the  sure-kill  temperature  of  the  skin  material,  then  the 
probability  of  survival  will  be  zero.  If  the  temperature 
is  below  the  sure-safe  temperature,  the  probability  of  sur¬ 
vival  will  be  one. 

In  current  multiburst  scenarios,  the  thermal  effect 
due  to  each  burst  is  treated  as  an  independent  event,  mean¬ 
ing  that  the  temperature  rise  of  the  target  skin  is  con¬ 
sidered  to  be  the  result  of  thermal  radiation  emitted  from 
only  one  fireball.  For  a  series  of  bursts,  the  temperature 
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rise  would  be  monitored  individually,  resulting  in  a  proba¬ 
bility  of  survival  for  each  burst.  The  probabilities  would 
then  be  combined  to  arrive  at  a  final  probability  of  sur¬ 
vival  for  the  target.  This  method  of  treating  bursts  as 
independent  events  will  be  called  the  noncumulativa  case. 
The  cumulative  case  is  the  case  of  considering  the  tempera¬ 
ture  rise  of  the  target  skin  to  be  the  result  of  receiving 
thermal  radiation  from  more  than  one  fireball  at  the  same 
v.ime . 

Problem  and  Scope 

The  purpose  of  this  thesis  project  was  to  develop 
a  computer  program  to  calculate  the  skin  temperature  rise 
of  a  missile  subjected  to  cumulative  thermal  effects  from  a 
series  of  bursts  occurring  during  the  missile's  flight. 

The  probability  of  survival  for  this  case  was  then  compared 
to  the  probability  calculated  for  the  noncumulativa  case  to 
demonstrate  the  need  for  considering  cumulative  effects. 
Furthermore,  since  blast  effects  are  generally  considered 
more  lethal  than  thermal  effects  in  the  noncumulative  case, 
the  cumulative  results  were  also  compared  to  results  from 
noncumulative  blast  effects  to  determine  if  thermal  radia¬ 
tion  becomes  the  kill  mechanism.  During  this  research, 
no  effort  was  made  to  vary  the  threat  conditions  or  to 
optimize  the  probability  of  survival  for  a  missile. 

Instead,  the  intent  was  to  show  that  treating  each  burst 
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as  an  independent  event  gives  results  that  severely  over- 
estimate  a  missile's  probability  of  survival. 

The  burst-target  missile  system  studied  was  the 
Peacekeeper  close-spaced  basing  (CSB)  system,  more  commonly 
known  as  Dense  Pack.  Although  this  system  is  no  longer 
being  considered  for  national  defense,  it  is  useful  for 
showing  the  differences  between  the  cumulative  and  noncumula- 
tive  cases.  The  configuration  of  the  missile  field  is  shown 
in  figure  1.  The  field  is  subjected  to  a  walk  attaak  start¬ 
ing  at  silo  #1  and  proceeding  every  two  seconds  to  silos 
#2,  #3,  #4,  and  so  on.  The  basic  scenario  studied  was  that 
of  a  missile  launching  at  the  same  time  as  silo  #1  is  hit. 

Assumptions  and  Limitations 

The  following  assumptions  or  limitations  were  made 
to  simplify  the  analysis i 

1.  Cratering  and  blast  effects  produoed  by  the 
burst  are  not  modeled. 

2.  The  fireball  is  considered  to  be  an  isotropio 
point  source  of  thermal  radiation.  Furthermore,  the  point 
source  is  assumed  to  be  centered  at  the  top  of  the  rising 
dust  cloud  created  by  the  fireball. 

3.  The  amount  of  scatter  and  absorption  of  thermal 
radiation  in  the  air  is  quantified  by  using  the  trans¬ 
mittance  t,  which  is  the  fraction  of  direct  and  scattered 
radiation  transmitted  from  burst  to  target.  An  expression 
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for  t  as  a  function  of  slant  range  will  be  developed  from 
data  presented  for  a  cloudless  atmosphere  with  a  visibility 
of  ~20  km  (Glasstone  and  Dolan,  1977:318).  This  expression 
for  t  does  not  account  for  the  decrease  in  transmitted 
radiation  due  to  the  dust  that  has  been  lifted  into  the 
air. 

4.  The  thermal  conductivity  and  specific  heat  of 
the  missile  skin  material  are  assumed  to  be  independent  of 
temperature  and  therefore  constant  during  the  heating  of 
the  skin.  This  is  a  reasonable  assumption  for  temperatures 
between  619  °K  and  809  °K  (Touloukian  et  al. ,  1970 <  1 

(Vol  1)  and  1  (Vol  4)],  the  sure-safe  and  sure-kill  inten¬ 
sities  chosen  for  an  aluminum  missile  skin.  Although  the 
assumption  is  not  true  for  temperatures  above  the  melting 
point  of  the  missile  skin,  the  change  in  material  proper¬ 
ties  is  not  important  compared  to  the  fact  that  melting 
constitutes  a  1001  failure  of  the  missile. 

5.  The  missile  skin  surface  is  modeled  as  a  flat 
plate  rather  than  as  its  actual  shape  of  a  cylinder.  This 
assumption  is  reasonable  for  an  energy  balance  made  over  a 
square  meter  of  skin  surface.  For  a  Peacekeeper  missile 
with  a  diameter  of  2.3  m  (Young,  198 4 1 168 )  ,  the  difference 
between  a  flat  surface  and  a  curved  surface  is  ~3%,  as 
shown  in  figure  2. 

Also,  the  slant  range  from  burst  to  missile  is  not 
calculated  for  a  specific  point  along  the  missile's  length. 
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Again,  for  a  Paacakaapar  miaaila  with  a  langth  of  *21  m 
(Young,  1984* 168) ,  tha  diffaranca  in  alant  rangaa  calcu- 
latad  from  aithar  and  of  tha  miaaila  would  not  ba  aignifi- 
cant  for  tha  diatanoaa  involvad. 

6.  Tha  miaaila1 a  poaition  and  valooity  at  any 
tima  ara  givan  in  figura  3.  Tha  miaaila  haa  no  valooity 
componant  in  tha  y  diraotion  and  tharafora  fliaa  atraight 
north,  aa  indicatad  in  figura  1. 

7.  Tha  walk  attack  ia  aaaumad  to  occur  with  par- 
fact  timing  and  at  a  parfaot  haight  of  burst  of  0  m  (i.a., 
a  aurfaca  burat) .  Thua,  tha  only  aiming  arror  oonaidarad 

ia  tha  RV'a  horizontal  distance  from  daaignatad  ground  zaro, 


where  designated  ground  zero  is  the  location  of  the  tar¬ 
geted  silo. 

8.  The  probability  of  damage  is  a  function  of 
intensity  rather  than  range  and  is  described  using  a  cumula¬ 
tive  log-normal  distribution  function. 

9.  The  analysis  considers  only  bursts  that  occur 
on  or  after  a  missile's  launch  time.  This  assumption 
simplifies  the  computer  simulation  of  the  problem. 

Approach 

An  expression  for  the  rate  of  heating  at  the  sur¬ 
face  of  the  missile  skin  was  derived  from  an  energy  balance 
over  a  square  mater  of  the  flat  surface.  This  balance 
included  absorption,  convection,  and  re-radiation,  but  it 
was  later  determined  that  re-radiation  was  not  significant 
for  the  temperatures  of  interest.  The  energy  balance  was 
than  solved  for  the  missile  skin  temperature  using  a  simpli¬ 
fication  called  the  thin  skin  approximation.  The  calcula¬ 
tion  of  the  maximum  rikiri  temperature  for  a  series  of  bursts 
involved  an  iterative  process  that  modeled  the  cumulative 
offeot  of  more  than  one  burst  by  combining  the  amount  of 
thermal  radiation  emitted  by  each  burst  at  a  particular 
time.  The  maximum  skin  temperature  was  then  used  to  calcu¬ 
late  the  probability  of  survival  for  the  missile. 
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Presentation 


The  derivation  of  the  equation  for  the  missile  skin 
temperature,  a  discussion  of  how  the  probabilities  were 
calculated,  and  an  explanation  of  how  this  information  was 
used  to  determine  the  probability  of  survival  in  both  the 
cumulative  and  noncumulative  cases  is  presented  in  Chap¬ 
ter  II.  Chapter  III  suipmarizes  the  conditions  and  param¬ 
eters  used  in  this  study  and  explains  the  reasoning  behind 
the  choice  of  the  parameters.  The  results  of  the  study  are 
discussed  in  Chapter  IV  and  conclusions  and  recommendations 
are  presented  in  Chapter  V. 
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II.  Theory 


This  ohapter  contains  the  theory  used  to  determine 
the  probability  of  survival  of  a  missile  subjected  to 
thermal  radiation  in  both  the  cumulative  and  noncumulative 
cases.  First,  the  derivation  of  the  equation  for  the  maxi¬ 
mum  skin  temperature  is  presented,  followed  by  an  explana¬ 
tion  of  how  this  equation  is  used  for  a  single  burst.  Next, 
the  probability  of  damage  function  and  aiming  error  are 
introduced.  This  information  is  then  used  to  explain  how 
the  probability  of  survival  is  calculated  for  the  noncumula- 
tive  and  cumulative  case. 

Derivation  of  the  Equation  for 
Maximum  sic  in  Temperature 

The  missile's  probability  of  survival  for  the 
thermal  threat  is  calculated  using  the  maximum  missile  skin 
temperature  reached  during  a  burst  scenario.  An  equation 
for  this  temperature  is  derived  using  an  energy  balance 
over  a  unit  are*  of  missile  skin  surface  and  a  simplifica¬ 
tion  known  as  the  thin  skin  approximation.  The  following 
presentation  of  the  derivation  is  taken  from  McKee  (McKee, 
1984:1-8)  . 

The  cylindrical  shape  of  the  missile  is  not  con¬ 
sidered  for  this  derivation.  Instead,  the  missile  skin  is 
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modeled  as  a  finite  slab  of  thickness  d,  with  no  heat  lost 
from  the  back  wall  of  the  skin  (x  -  d)  and  no  heating  from 
air  friction.  The  energy  balance  for  a  unit  area  of  missile 
skin  surface  (x  ■  0)  in  terms  of  fluence  is: 

Fheat  “  F absorbed  ~  ^convection  ”  Fradiation  (2.1) 

where 

2 

fcV  .  ■  Total  amount  of  energy /m  available  to 
n  heat  the  missile  skin  at  x  -  0 

'  absorbed  -  Amount  of  incident  thermal  radiation 
that  is  absorbed  by  the  missile  akin 

_  ■  Amount  of  absorbed  radiation  that  is 
convection  lost  to  convective  cooling  by  air 

P  nj.u..  ■  Amount  of  absorbed  radiation  that  is 
raaiat  re-radiated  as  1R  black-body  radiation 

to  the  environment 

Since  this  equation  must  hold  for  all  times,  the 
time  derivative  may  be  obtained  to  yield  the  following 
differential  equation: 

Fheat  “  Fabaorbed  "  1  convection  ~  Fradiation  (2.2) 


where  the  dots  indicate  the  time  derivative  or  a  rate. 

The  total  amount  of  energy  available  for  heating, 
Fheat'  k®  conducted  from  the  skin  surface  through  the 

skin  thickness.  If  the  flow  of  heat  is  assumed  to  be  one¬ 
dimensional  and  the  properties  of  the  medium  are  constant, 
then  the  conduction  process  is  described  by  the  heat 
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transfer  equation  (Holman,  1976sl02): 


3 


1  3T (x,  t) 
K  3t 


(2.3) 


where 

K  ■  thermal  diffusivity  (m2/sec)  ■  -~ 

«pP 

k  ■  thermal  conductivity  (J/m~s-°K) 
c  «  specific  heat  capacity  (J/kg-°K) 

hr 

2 

p  »  density  (kg/m  ) 

T(x,t)  ■  skin  temperature  (°K)  at  depth  x  and  time  t 

This  equation  can  be  used  to  estimate  the  thermal  diffu¬ 
sion  time  tdiff  required  for  heat  to  be  conducted  through 
the  skin  thickness  d  (see  Appendix  A).  The  result  1st 

^diff  ™  ^  ^  (2.4) 


If  tdi£f  is  much  less  than  the  time  scale  over  which  the 

thermal  pulse  occurs,  the  missile  skin  will  essentially  be 

at  a  uniform  temperature  throughout  x  at  any  given  time  t. 

The  thermal  pulse  time  scale  is  t  ,  the  time  of  the  second 

max 

thermal  maximum.  For  air  bursts  below  4572  m,  the  expres¬ 
sion  for  t  is  (Glasstone  and  Dolan,  1977s310) i 
max 

t  -  0.0417  Y*44  (2.5) 

max 


where  Y  is  the  yield  in  kilotons.  This  expression  will 
also  be  used  for  surface  bursts.  Thus,  for  uniform  heating 
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to  occur: 


*di*f  "  W  or  d  "  <K,W 


.5 


(2.6) 


If  d  fulfills  the  above  requirement,  then  the  thin  skin 
approximation  of  uniform  temperature  throughout  the  skin 
thickness  is  valid.  As  an  example,  tmax  is  "1.18  sea  for 
a  2  NT  burst  and  the  thermal  diffusivity  of  aluminum  is 
1  cm2/sec,  so  that  (K*t_,_w)*5  is  "1.0  cm.  Therefore,  a 
reasonable  value  for  d  would  be  ten  times  less  than  this 
value,  or  0.1  cm. 

An  expression  for  Fh#at  in  equation  (2.2)  oan  be 
derived  from  equation  (2.3),  the  thin  skin  approximation, 
and  Fourier's  law  of  heat  conduction  (Holman,  1976 1 2) s 


F(x,t) 


gJJfeft) 

*  ax 


(2.7) 


where  F  is  the  rate  of  heating.  Differentiating  (2.7) 
with  respect  to  x  and  substituting  into  (2.3)  yields: 


ax 


v3^ 


(2.8) 


If  the  thin  skin  approximation  applies,  then  T(x,t)  becomes 
a  function  of  time  only  and  the  temperature  at  x  -  0  equals 
the  temperature  at  x  ■  d  equals  T(t).  The  above  equation 
can  then  be  integrated  over  the  skin  thickness: 


-C  ^  -  -//  v 


(2.9) 
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to  obtains 


F(0,t)  -  F(d,t)  -  oppd  (^5^-) 


(2.10) 


By  the  definition  of  and  the  assumption  that  no  heat 

is  lost  at  x  *  d,  the  final  result  iss 


F. 


heat 


* 


(2.11) 


where  a  ■  c  pd. 

P 


The  amount  of  incident  radiation  absorbed  ist 


'absorbed  " a  F incident 


(2.12) 


where 


a  -  absorptivity  of  missile  skin 


radiant  exposure  from  a  burst  (Glasstona  and 
incident  Dolan<  1977,316) 


,12. 


.  CP  ,  tf.y^.isexio-jT 

4^ (SR) 2 


(2.13) 


tf  -  thermal  fraction  or  effective  thermal  par¬ 
titions  fraction  of  bomb  yield  appearing  in 
the  form  of  thermal  radiation  ■  . 18  for  a 
surface  burst  (Glasstone  and  Dolan,  1977s 319) 


y  ■  bomb  yield  (kt) 
t  ■  transmittance  of  atmosphere 

SR  -  slant  range  from  thermal  radiation  point 
source  to  missile  (m) 

CF  ■  a  correction  factors  sin(^) ,  where  ^  ■  angle 
between  missile  skin  surface  and  slant  range 
vector 
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The  energy  absorption  rate  ist 


absorbed 


incident 

dt 


(2.14) 


The  convective  cooling  term  can  be  written  immedi¬ 
ately  as  a  rate  as  follows  (Holman,  1976tl2)t 


^conviction  “  h[T<t)  '  Vr*4'1 


(2.15) 


where 


-  local  convs 
(J/m2-«-°K) 


convective  heat  transfer  coefficient 


Tair  (t» 


temperature  of  ambient  air  at  missile  altitude 
(°K) 


The  variable  h  depends  on  the  velocity  and  temperature  of 
the  air  flowing  along  the  missile  and  is  calculated  for  a 
specific  point  on  the  missile  skin.  Appendix  E  describes 
how  h  is  calculated  for  a  flat  plate  heated  to  a  uniform 
temperature  over  its  entire  length. 

The  expression  for  the  rate  of  black-body  radiation 


is  (Holman,  197 6 s 13 ) t 


radiation 


o[T(t)  4  -  T.ir  <t)4] 


(2.16) 


where  a  is  the  Stefan-Boltzmann  constant  -  5.6696X10 
J/m2-s-8K4.  The  radiation  term  will  be  insignificant  for 
skin  temperatures  below  1000  °K,  which  is  the  case  for  the 


•  .  If  i*  iw, 

•  .  4  .  «*<•  #'* . V  «. •  w.a '•» 


Is 


scenarios  considered  in  this  study.  Thus,  radiation  is 
neglected  in  the  analysis. 

Equation  (2.2)  can  now  be  written  asi 

dF 

a  m  ot"-in’dtg,nt  ”  h[T(t)  “  Tair{t)]  U*17) 

This  differential  equation  is  solved  using  the  method  of 
finite  differences.  Choosing  t  v  an  a  finite  time  step, 

IuCmw 

the  following  replacements  are  madet 

Ml  +  hlh. 

^  Snax 


^incident 
dt  “ 


-AS- 

Snax 


T  (t) 


where  the  time  steps  are  numbered  as  j  *  1,2,3,...,N  andi 

T,  “  temperature  at  the  end  of  the  jth  time  step 
2  (°K) 

T,  ■  temperature  at  the  beginning  of  the  jth  time 
A  step  ( °K) 

AQ  «  total  thermal  fluence  that  is  incident  on  the 
missile  during  the  jth  time  step  (J/m2) 


Specifically* 


-  aCTj  •  ^incident 
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where  pinci(jent  i-  <*«fined  in  equation  (2.13)  and  ACT^  ia 
the  difference  between  the  fraction  of  thermal  energy 
emitted  up  to  the  beginning  and  up  to  the  end  of  Lhe  jth 
time  step.  The  fraction  of  thermal  energy  emitted  at 
time  t,  CT j i  is  plotted  versus  normalized  time  a 

the  right  hand  curve  in  figure  4. 

After  making  the  above  replacements  and  solving 
for  T,,  equation  (2.17)  becomes » 


Thus,  the  missile  skin  temperature  at  the  end  of  time 
step  j  can  be  found  using  equation  (2.19).  This  equation 
is  next  used  to  calculate  the  maximum  akin  temperature 
reached  for  a  single  burst. 


MUldM-.MJlil’l 


the  Maximum  Skin 


The  following  definitions  are  needed  to  describe 
how  the  maximum  akin  temperature  is  calculated  for  a  single 
burst i 


t  •  time  of  missile  skin  exposure  to  thermal 

energy,  t  is  also  equal  to  the  time  of  thermal 
energy  emission  by  a  burst  because  radiation 
travels  at  the  speed  of  light. 

t  -  time  into  missile  flight 
111 

It  ■  launch  time  of  the  missile.  At  It,  tm  -  0 
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Fig-  4.  Fraction  of  Thermal  Energy  vs.  Normalized  T 
(Glasstone  and  Dolan,  1977:311) 


4> 


tb  -  time  of  bunt.  At  tb,  t  *  0 

tn  ■  tima  into  miaaila  flight  whan  burat  ocoura 
0  tQ  -  tb  It. 

Thua,,  two  indapandant  tima  lima  a  ira  occurring*  ona  dafining 
tha  miaaila  flight  and  ona  dafining  tha  tharmal  radiation 
amittad  by  a  burat.  Figura  5  ahowk?  tha  ralation  batwaan 
thaaa  tima  linaa  and  how  j  ia  dafinad  for  aaoh  tima  atap 
aftar  a  burat. 


t  •  tint  of 
axposura  to 
tharmal  radi¬ 
ation 

0 

1 

2 

3 

j 

t *  tima  into 

♦ 

1 

H— 

B  miaaila 

t 

1 

t- 

1 

1 

1 

rt 

o 

-l-o 

ltaaa 

2tmax 

3t 

MX 

flight 

It*  launch  tima 

tb*  tima  of  burat 

It 

tb 

t_+2t 

0  max 

tft+3t  "  B 

0  aax 

tQ*  tb-lt 
j  ■  tima  atap  # 

Fig.  S.  , singla  Burat  Tima  Linaa 


Tha  prooadura  for  datarmining  tha  maximum  akin 
tamparatura  ia  an  itarativa  prooaaa  that  bagina  by  calcu- 
lating  at  tha  and  of  tha  firat  tima  atap  (j  ■  1) .  To 
do  thia  calculation,  tha  variablaa  T^,  AO,  Tair,  and  h 
muat  ba  known.  For  tha  firat  tima  atap,  T^  ia  aaaumad  to 
ba  aqual  to  Tt^r  at  tha  miaaila  altituda  at  tha  tima  of 
burat.  Tha  actual  valua  of  would  ba  highar  than  Tair 
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due  to  aerodynamic  heating,  but  the  aasumption  is  good  for 
a  elowly  moving  miaaile.  AQ,  Taij.,  and  h  are  dependent  on 
miaaile  and/or  burat  oharaoteriatica  that  change  with  time. 
The  miaaile  oharaoteriatica  are  velocity,  altitude,  down- 
range-diatanoe,  and  the  flight  path  angle  9,  the  angle 
between  the  velocity  vector  and  the  horizontal  plane.  The 
burat  oharaoteriatica  are  CTj  and  the  height  of  the  fire¬ 
ball  Hfbj .  If  the  fireball  ia  aa Burned  to  be  at  the  top  of 
the  duat  cloud,  the  following  equation  can  be  uaed  to  oalau- 
late  Hfbj  (McGahan  et  al.,  1971t4Q)i 

ttfbj  -  21,«4Q.8(W*177)  tl-U-t/t^2]  (2.20) 

where  W  *  yield  in  megatona  and  tg  -  aloud  atabilization 
time  -  240  aec  (Bridgman,  1984). 

Miaaile  char ac ter iatioa  are  calculated  at  an 
average  time  tm  -  tQ  +  (j-.S)*tmax,  the  midpoint  of  each 
time  atep.  The  altitude  and  velocity  are  then  uaed  to 
calculate  Tftir  and  h.  Thia  value  for  h  repreaenta  the 
average  amount  of  convective  cooling  that  occurred  during 
the  jth  time  atep.  To  calculate  the  alant  range  (SR) ,  all 
miaaile  oharaoteriatica  are  needed  and  the  height  of  the 
fireball  muat  be  known  at  time  t  ■  (j-.5)*t  .  Once  SR 

ia  known,  the  variablea  CF  and  t  can  be  determined.  The 
quantity  ACT j  ia  calculated  for  the  time  between  j  and  j-l, 
and  ia  uaed  with  SR,  CF,  and  t  to  determine  AQ.  Detaila  on 
how  to  calculate  miaaile  oharaoteriatica,  burat 
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characteristic*,  h,  ta^r,  SR,  CF,  and  ?  are  found  in 
Appendices  B,  C,  D,  E,  and  F. 

Once  the  necessary  variables  are  known,  t2  can  be 
aaloulated  for  time  step  j.  For  the  next  time  step, 
is  set  equal  to  T2  and  new  variables  are  determined  to 
calculate  a  new  T2  for  that  time  step.  This  iteration 
process  is  continued  until! 

1.  t2  is  found  to  be  less  than  T^.  This  indi¬ 
cates  that  the  amount  of  heat  removed  by  convective  cooling 
is  now  greater  than  the  amount  absorbed  and  therefore  the 
skin  temperature  will  continue  to  decrease.  The  maximum 
skin  temperature  reached  is  the  current  value  of  T^t  or 

2.  The  number  of  time  steps  j  equals  10.  The 
process  is  stopped  here  because  801  of  the  thermal  radia¬ 
tion  from  the  burst  has  been  emitted,  and  the  value  of 
ACTj  would  be  negligible  for  the  remaining  time  steps  (see 
figure  4).  The  maximum  skin  temperature  reached  is  the 
value  of  T2  at  j  -  10. 

The  iteration  process  to  calculate  the  maximum  skin  tempera¬ 
ture  for  a  single  burst  is  summarised  in  the  algorithm 
prmrnted  in  table  I.  The  maximum  skin  temperature  is  then 
used  to  determine  the  probability  of  damage. 

Curvivabillty  and  Aiming  Error 

For  th.'.c  analysis,  the  probability  of  damage  is 
based  on  the  >kre  nafe  and  sure  kill  intensities  laB  and  I 
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TABLE  I 


ALGORITHM  FOR  FINDING  MaxT  FOR  A  SINGLE  BURST 


1.  Sat  j  ■  1,  Tj_ 


Tair  at  ^ 


2. 


Find  miaaile  velocity, 
distanca,  and  e  at  < 
(Saa  Appendix  B) 


altitude,  down-range* 

•  t*  +  (j-.5)*t 
0  'J  1  max 


Find  fireball  height  at  t  ■  (j-.5)*t 
(See  equation  (2.20))  max 


Find  ACTj  for  the  time  j-1  to  j 
(See  Appendix  F) 


3.  Calculate  h,  Tft^_ ,  SR,  CF,  and  t 
(See  Appendioea  C,  D,  and  E) 


4.  Calculate  AQ  uaing  equation  (2.18) 


5.  Calculate  T2  uaing  equation  (2.19) 

6.  If  Tj  <  T]_,  then  go  to  atep  8 
Otherwiae,  aet  T.^  ■  T2 

7.  If  j  <  10,  then  j  ■  j  +  1  and  go  to  atep  2 
Otherwiae,  go  to  atep  8 

8.  Maximum  akin  temperature  MaxT  ■  T^ 


where  intensity  refers  to  temperature.  For  aluminum, 

I  and  I  ,,  are  chosen  to  be  619  #K  and  809  °K  respee- 
tively  (Bridgman,  1984) .  The  probability  of  damage  is 
calculated  using  the  cumulative  log-normal  distribution 
function! 


Pd(D 


1 

/57s 'i 


In  1  - 

V 


a 


.  2 

)  ]dl 


(2.21) 


where  the  intensity  I  refers  to  the  maximum  skin  tempera¬ 
ture,  MaxT.  The  parameters  a'  and  S'  are  calculated  using 
I  and  IbJc.  By  definition,  if  Pd(I>g)  "0.98  and  Pd(IgJc) 

■  0.02,  then  from  equation  (2.21) t 

■  2  l"  <2-221 


3' 


T7OT 


In  («■*) 

ss 


(2.23) 


A  more  detailed  discussion  of  P^d)  ,  a',  and  3'  is  pre¬ 
sented  in  Appendix  G. 

The  probability  of  damage  calculated  using  the 
temperature  determined  from  the  procedure  in  table  I  is 
the  probability  of  damage  from  an  RV  that  was  assumed  to 
hit  at  its  designated  ground  zero,  or  at  the  center  of  the 
targeted  silo.  However,  in  reality  each  RV  will  have  an 
aiming  error  that  could  cause  it  to  land  at  a  point  other 
than  its  designated  ground  zero.  Since  the  position  of 
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the  burst  must  be  known  to  calculate  MaxT,  the  aiming  error- 
needs  to  be  quantified.  This  is  done  by  dividing  the  area 
around  a  silo  into  discrete  cells  whose  dimensions  are 
such  that  each  cell  has  an  equal  probability  of  being  hit. 
Although  the  dimensions  of  these  cells  are  fixed,  the  posi¬ 
tion  of  the  cells  around  the  designated  ground  zero  is  not 
unique.  Figure  6  shows  the  ten  cell  configuration  used  in 
this  analysis.  It  also  shows  <p^>  and  <6^>,  the  two  param¬ 
eters  that  locate  the  call  centroids  such  that  each  cell  i 
is  of  equal  probability.  The  values  of  <p^>  depend  on  the 
circular  error  probable,  or  CEP,  which  is  defined  as  the 
radius  inside  of  which  50%  of  the  incoming  RVs  will  hit. 

The  method  of  calculating  <p^>  and  <6^>  is  described  in 
Appendix  G. 

The  point  defined  by  <P1>  and  <0^>  for  cell  i  is 
used  as  ground  zero  for  calculating  MaxT.  Thus,  for  a 
•ingle  burst,  tan  different  burst-missile  configurations 
are  possible  and  each  configuration  will  have  an  associ¬ 
ated  probability  of  damage.  The  total  probability  of 
damage  for  N  cells  ist 

1  N 

P-  -  w  2  P-  (MaxT)  (2.24) 

d  N  i-1  di 

and  the  probability  of  survival  is  then: 

P8  -  1  -  Pd  (2.25) 
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Fig.  6.  Ten-cell  CEP 

The  calculation  of  P#  for  a  aeries  of  bursts  depends  on 
whether  the  noncumulative  or  cumulative  case  is  being  con¬ 
sidered. 

Calculating  the  Probability  of 
Survival  £ or  a  Noncumulative 
Burst  Scenario 

For  the  noncumulative  case,  each  burst  is  treated 
as  an  independent  event.  Thus,  the  temperature  rise  caused 
by  one  burst  is  not  affected  by  the  additional  thermal  radi¬ 
ation  from  a  subsequent  burst.  Each  independent  event  k  has 
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an  associated  probability  of  survival  Pg  calculated  from 
equation  (2.25).  The  total  probability  of  survival  for  all 
bursts  considered  is  then: 


Total  m 


total  # 
of  bursts 


(2.26) 


calculating  the  Probability  of 
Survival  for  a  Cumulative 


Survival  tor  a  Cumulative 
Burst  Scenario  ” 


For  the  cumulative  case,  each  burst  is  not  treated 
as  an  independent  event.  Rather,  the  amount  of  thermal 
energy  emitted  by  each  burst  will  affect  T2  through  the  term 
AQ  in  equation  (2.13).  Equation  (2.18)  can  be  rewritten  as 
an  expression  for  AQ  for  a  burst  k: 


,12  / CFkrk  ACT.k 


l  UV*  J 

(SRk) 2 


(2.27) 


where  y  is  the  same  for  each  burst  and  the  terms  in  paren- 
theses  are  burst  dependent.  SR  ,  CF  ,  and  t  are  calculated 
using  the  position  of  burst  k,  and  ACTj  is  calculated  based 
on  the  time  of  burst  k. 

Figure  7  shows  the  timing  for  a  series  of  four 
2  megaton  bursts.  The  missile  is  launched  at  It,  the  first 
burst  occurs  at  tb^,  and  all  other  bursts  occur  at  two 
second  intervals  after  tb^.  All  calculations  for  bursts 
2,  3,  and  4  are  referenced  to  the  time  steps  of  burst  1. 


'  -rime  am  rtii  vns  inr 1  “r  i  ret 


Multiburst  Tine  Lines 


For  2  MT  weapons,  t  -1.18  seconds,  so  the  time  steps 

max 

are  shown  to  be  shorter  than  the  interval  between  bursts. 

Because  thermal  energy  is  additive,  the  term  Aq  in 

equation  (2.19)  will  be  a  sum  of  the  thermal  energy  emitted 

by  each  burst  that  has  been  detonated  by  time  t.  For 

example,  at  j  -  3,  the  missile  will  be  exposed  to  the 

thermal  energy  from  two  bursts.  Aq1  will  be  calculated  as 

described  in  table  I,  using  missile  characteristics  and 

Hfbj1  calculated  at  the  midpoint  of  the  time  step,  and 

Act.1  calculated  for  the  time  between  j  -  2  and  j  -  3.  For 
3  2 

the  calculation  of  A q‘  at  j  -  3,  the  missile  characteris¬ 
tics  will  be  the  same  as  those  for  Aq1  since  they  are  based 
on  missile  launch  time  and  therefore  do  not  depend  on  burst 
time.  However,  the  slant  range  will  be  different  because 

of  the  different  position  and  fireball  height  of  burst  2. 

2 

The  quantity  ACT^  will  also  be  different  because  burst  2 

has  not  emitted  as  much  radiation  as  burst  1  for  the  time 

between  t/t  -  2  and  t/tm,„  -  3. 
max  max 

For  a  burst  k  occurring  at  time  tb^  (whore  k  >  1) , 

A  CT  .k  will  be  calculated  at  time  t"  -  (tb,  +  j*t_,_„ 

3  l  max 

-  tb.)/tm-v.  For  the  time  steps  where  the  emission  of  a 
x  max 

burst  k  (k  >  1)  does  not  last  the  full  time  step,  as  for 
j  »  2,  4,  and  6,  the  calculation  of  Aq'v  is  slightly  altered 
by  finding  the  missile  characteristics  at  a  time  that  is 
the  midpoint  of  the  first  "shortened"  time  step  of  burst  k. 
Once  all  the  AQ  ■  s  are  calculated,  they  are  added  to  get 
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the  total  AQ,  and  this  term  Is  then  used  In  equation 
(2.19)  . 

Total 

The  calculation  of  P  also  changes  for  the 

JB 

cumulative  case.  Instead  of  a  series  of  four  10-cell  con¬ 
figurations,  now  each  combination  of  configurations  must 
be  considered.  For  k  bursts,  each  with  a  10-cell  CEP,  a 
total  of  10  configurations  are  possible  and  therefore 

V 

10  Pd(MaxT) 's  must  be  calculated.  The  total  probability 
of  survival  for  k  bursts  is  then* 

Total  m  l  m  l  L  p  (MaxT)  (2.28) 

*  10*  i-1  ai 

The  equation  for  the  maximum  temperature  and  the 
information  on  how  to  calculate  the  probability  of  survival 
oan  now  be  used  to  develop  a  computer  program  that  deter¬ 
mines  P0  for  a  series  of  bursts.  However,  the  final  ver¬ 
sion  of  this  program  depends  on  parameters  that  are  dis¬ 
cussed  in  the  next  chapter. 
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III.  Problem  Parameters 

The  solution  to  the  problem  of  determining  the 
maximum  akin  temperature  and  the  probability  of  survival  for 
a  missile  required  a  number  of  specific  burst  and  missile 
conditions  and  also  the  choice  of  certain  limiting  param¬ 
eters  .  These  conditions  and  parameters  are  summarized  in 
table  II.  The  remainder  of  this  chapter  explains  the 
reasons  for  choosing  the  parameters  listed  in  table  II  and 
also  describes  how  the  calculation  of  the  maximum  skin 
temperature  was  adapted  to  solution  on  a  computer. 

Choosing  the  Maximum  Number 
oT  Bursts 

The  maximum  number  of  bursts  affecting  the  missile , 
maxb,  is  an  important  quantity  because  it  determines  the 
number  of  burst-missile  configurations  that  need  to  be  con¬ 
sidered  in  order  to  calculate  the  probability  of  survival 
for  the  missile.  For  each  burst  with  a  10-cell  CEP,  a  total 
of  10maxb  configurations  are  possible  and  therefore  iom<lxb 
maximum  skin  temperatures  must  be  calculated. 

Because  the  missile  is  accelerating  and  rising  away 
from  the  fireballs ,  it  is  possible  that  an  upper  limit  to 
the  number  of  bursts  exists.  Any  bursts  occurring  after 
this  limit  would  not  raise  the  skin  temperature  further. 


TABLE  II 


CONDITIONS  AND  PARAMETERS  CHOSEN  FOR  STUDY 


Svatam i 

cToee-Spaced  Baaing,  or  Danaa  Pack  (aaa  figura  1) 
Thraat  Conditional 

Walk  attack  atarting  at  ailo  #1  and  continuing  avary 
2  aaoonda  on  auooaaaiva  ailos 

Waapon  Yield!  2  MT 

Haight  of  Buret t  0  m 

For  aurfaoa  burata,  tf  -  .18 
(Glaaatone  and  Dolan,  1977 t 319) 

Mlaalla  Conditional 

Mi'aai  la  via  loci  ty ,  altitude,  dovm-range-dietnnce,  and 
flight  path  angle  ahown  in  figure  3  aa  a  function 
of  time 

skin  material i  Aluminum 
K  *  .0001  m2/a 
p  -  2700  kg/m3 
On  «  900  J/kg-°K 
h  -  .50 
Igi  -  619  #K 
m  809  #K 

Skin  thioknaaai  d  -  .001  m 


Probability  Conditional 

RV  aiming  error  only  for  ground  aero  10-call  CEP 
Probability  of  damage  baaed  on  intanaity  and  calcu¬ 
lated  uaing  the  cumulative  log-normal  diatribution 
function 


Pd#daa) 

pd(I.k> 


.98 

.02 


Parameter a t 

Maximum  number  of  burata  conaideredi  maxb  «  4 

Maximum  number  of  time  a tape  needed!  11 

Heat  tranafer  coefficient  calculated  at  x  «  5.5  m 


This  uppar  limit  can  be  estimated  quickly  for  any  missile 
by  calculating  ?2  for  a  series  of  bursts  and  assuming  that 
each  RV  lands  at  its  designated  ground  zero  (i.e.,  at  the 
center  of  cell  #1  or  a  1-cell  CEP) .  Appendix  H  shows  that 
the  estimated  upper  limit  to  the  number  of  bursts  was  found 
to  be  eight,  requiring  the  maximum  akin  temperature  to  be 

o 

calculated  10  times.  However,  in  the  interest  of  reducing 
the  amount  of  computer  time,  the  maximum  number  of  allowed 
bursts  was  chosen  to  be  four.  Although  the  probability  of 
survival  for  a  four-burst  scenario  is  an  optimistic  result, 
the  scenario  adequately  demonstrates  the  effects  of  con¬ 
sidering  the  cumulative  versus  the  nonoumulative  case. 

Choosing  the  Maximum  Number 
of  Time  S topB~NeodeA 

For  the  cumulative  case,  computer  time  aan  be  saved 
if  the  number  of  time  steps  necessary  to  determine  the  maxi¬ 
mum  akin  temperature  is  known.  The  reason  for  this  will  be 
explained  later  in  this  chapter.  For  a  single  burst,  the 
time  considered  was  t  ■  10*tmax  (see  Chapter  II) .  However, 
for  multiple  bursts  it  is  possible  that  the  skin  tempera¬ 
ture  could  continue  to  rise  after  t  -  10*tmax,  even  though 
the  contribution  of  thermal  radiation  from  the  first  burst 
is  negligible  after  this  time. 

Figure  7  showed  the  duration  of  significant  thermal 
radiation  four  2  megaton  bursts  to  be  approximately  15*tmax 
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or  "17.8  seconds  long.  However,  Appendix  I  shows  that  for 
a  1-cell  CEP  situation,  the  maximum  skin  temperature  is 
reached  before  j  ■  10.  Therefore,  to  be  safe,  the  maximum 
number  of  time  steps  needed  was  chosen  to  be  eleven. 

Choosing  the  Position  for  the 
Local  Heat  Transfer  Coefficient 

The  last  decision  to  be  made  was  the  choice  of  xm, 
the  point  along  the  flat  plate  representing  the  missile  skin 
where  the  local  heat  transfer  coefficient  is  calculated. 
McKee  chose  x_  *  5.S  m,  the  location  of  the  third  stage 
joint  on  a  generic  missile  (McKee,  1984 s 7) .  Comparing 
values  of  h  for  various  choices  of  x]n  (sea  figure  E-l) 
shows  that  the  chosen  value  of  x^  gives  heat  transfer  coef- 
ficients  that  are  between  those  for  the  minimum  and  maximum 
amount  of  convective  cooling.  Thus,  the  location  of  xm 
-5.5m  represents  a  position  where  an  average  amount  of 
convective  cooling  occurs  along  the  missile  skin. 

Adapting  the  Problem  to 
a  Computer  Program 

As  explained  in  Chapter  II,  determining  the  proba¬ 
bility  of  survival  for  a  missile  exposed  to  k  bursts  in  the 

v 

cumulative  scenario  requires  calculating  10  maximum  skin 
temperatures,  one  calculation  for  each  geometric  configura¬ 
tion  due  to  the  10-cell  CEP  of  each  burst.  Of  the  vari¬ 
ables  used  to  calculate  Tj,  AQ,  or  h,  only  the  slant  range 
SR,  the  transmittance  x ,  and  the  term  CF  depend  on  the 
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buret  location.  Thus,  to  minimize  computer  time,  certain 
variables  can  be  calculated  once  and  stored  for  future  use. 

Specifically,  once  the  missile  launch  time  It,  the 
time  of  the  first  burst  tb^  and  the  total  number  of  neces¬ 
sary  time  steps  are  known,  the  midpoint  of  each  time  step  j 

can  be  calculated  as  t  «  t.  +  ( j-.  5)  *t  Missile  charaa- 

u  max 

terlstios,  the  ambient  air  temperature,  and  the  heat  trans¬ 
fer  coefficient  can  then  be  calculated  at  these  times  and 
stored  in  one-dimensional  arrays  whose  indices  refer  to  the 
particular  time  step  j.  Finally,  knowing  tho  total  number 
of  bursts  to  be  considered  and  the  time  of  eaoh  burst,  the 
terms  Act^  and  Hfb^  can  be  calculated  for  eaoh  time  step  j 
and  each  burst  k  and  stored  as  a  k  by  j  array. 

The  procedure  for  calculating  the  meximum  skin 
temperature  for  the  cumulative  case  can  now  be  presented 
using  information  from  tha  above  saotiona  and  Chapter  II. 
The  algorithm  is  shown  in  table  III.  As  mentioned  before, 
the  procedure  has  been  simplified  by  precalculating  and 
storing  certain  variables  for  use  daring  the  iterative 
proaess  of  calculating  T^. 

A  F0RTRAN7 7  computer  program  that  follows  tho 
algorithm  in  table  III  is  presented  in  Appendix  K.  This 
program  waa  written  on  a  Corona  PC  (MS-DOS)  and  was  also 
used  on  a  VAX  11-780  at  the  Air  Force  Institute  of  Tech¬ 
nology,  The  program  calculates  the  maximum  skin  tempera¬ 
ture  for  the  conditions  given  in  table  II  and  allown  the 
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TABLE  III 


ALGORITHM  FOR  FINDING  THE  MAXIMUM  SKIN  TEMPERATURE 
FOR  THE  CUMULATIVE  CASE 


1.  Knowing  It  and  tbi,  calculate  and  store  mil  'i la 

velocity,  altitude,  down-range-distance,  rnd 
6  at  tffl  *  tQ  +  (j-.5)*Tmax,  where  j  ■  1.  11. 

Caloulata  and  atora  T  and  h  at  tha  a.\m.'  timaa. 

Knowing  maxb  and  tho  timaa  of  tha  bursts,  calculate 
and  atora  Aq.  and  Hfbj  for  aach  burst  And  tint 
atap.  3  3 

2.  Sat  j  ■  1,  Tj  •  Tjiir  at  tm  -  to,  T2  to  ba  any 

number  greater  than  T. ,  ana  the  number  of  h  .rata 
nb  *  1.  1 

3.  If  Tj  >  T^  and  j  <  «  10  than* 

a.  Tha  currant  time  step  ia  j  ■  j  +  1>  AQ^  ■  0. 

b.  At  time  determine  if  another  burst  haa 

occurred.  If  so,  nb  *  nb  +  1  until  nb  *  4. 

c.  For  each  burst  k  that  has  occurred! 

1.  Calculate  SR,  GF,  and  t  using  tha  ground 
zero  determined  by  the  cell  f  of  burst  k, 
missile  characteristics  stored  at  time  j, 
and  Hfbj  for  burst  k.  If  burst  k  has  just 
occurred  within  tho  tima  atap,  missile  char¬ 
acteristics  must  ba  re-calculated  at  tha  mid¬ 
point  of  that  burst’s  first  shortened  time 
step. 

v 

2.  Calculate  AO  j  using  ACTj  for  burst  k, 
and  add  it  to J tha  currant  sum  of  . 

d.  Calculate  T2  knowing  h  and  tho  total  AO.,  for  time 

atop  j .  3 

e.  If  Tj  >  T:,  then  T2  «  T^ 

f„  Return  to  tho  condition  in  step  3.  If  either  test 

fails,  go  to  step  4. 

4.  MaxT  -  Tx. 
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1 

user  to  input  the  missile  number,  launch  time,  number  of 
bursts  to  be  considered  (<,  4)  ,  time  of  the  first  burst, 
number  of  cells  to  be  considered  (<  10) ,  and  the  option  of 
considering  rising  or  stationary  fireballs. 

Also  included  in  Appendix  K  is  a  BASIC  program  that 
calculates  the  maximum  skin  temperature  for  the  noncumula- 
tive  case.  This  program  follows  the  algorithm  given  in 
table  X  for  a  single  burst  and  does  not  need  to  preoalcu- 
late  or  store  any  values.  The  program  allows  the  user  to 
in$ut  the  same  quantities  as  the  program  written  for  the 
cumulative  case. 
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Chap tar  IV.  Result  a  and  Discussion 


Tha  main  objective  of  this  thesis  was  to  develop  a 
method  of  calculating  tha  maximum  skin  temperature  of  a 
missile  subjected  to  thermal  radiation  from  more  than  one 
fireball  at  a  time.  The  results  of  this  cumulative  case 
were  then  compared  to  the  results  of  the  noncumulativa  case, 
which  considered  the  effect  of  thermal  radiation  from  only 
one  fireball  at  a  time.  The  comparison  showed  that  missiles 
experienced  a  much  greater  temperature  rise  when  subjected 
to  radiation  from  four  fireballs  at  a  time  rather  than  when 
the  four  bursts  were  treated  as  independent  events.  Thus, 
the  more  realistic  case  of  cumulative  thermal  effects  indi¬ 
cates  that  targets  suffer  more  thermal  damage  in  a  multi¬ 
burst  scenario  than  previously  expected. 

The  results  of  the  cumulative  case  are  first  pre- 
‘  sented  for  a  1-cell  CEP;  i.e.,  for  a  perfect  hit  by  an  RV 
on  its  targeted  silo.  These  results  show  the  temperature 
rise  of  the  missile  skin.  Also  included  are  a  comparison 
of  considering  a  rising  fireball  at  the  top  of  a  dust  cloud 
to  a  stationary  fireball  on  the  ground,  and  the  effect  of 
decreasing  the  time  step  of  the  calculation  to  t_„v/2. 

Next,  the  cumulative  case  is  compared  to  the  noncumulative 
case  of  a  missile  subjected  to  four  bursts  with  a  1-cell 
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CEP.  This  comparison  shows  the  marked  difference  in  the 
maximum  skin  temperature  reached  in  both  cases.  The 
scenarios  are  then  extended  to  the  full  10-cell  CEP  for 
four  bursts  to  determine  which  silos  are  in  sure-kill  or 
sure-safe  territory.  The  10-oell  CEP  cumulative  thermal 
results  are  also  compared  to  noncumulative  blast  results 
to  determine  which  effect  is  more  lethal.  Finally#  two 
topics  related  to  the  cumulative  results  are  discussedt 
the  effects  of  dust  shielding  and  a  comparison  of  1-cell 
CEP  to  10-cell  CEP  results. 

Temperature  Rise  for  tho 
Cumulative  Case  (1 -cell  CEP) 

Figure  8  shows  the  temperature  rise  of  a  missile 
subjected  to  thermal  radiation  from  a  series  of  four  bursts. 
The  temperature  plotted#  T2#  is  the  temperature  at  the  and 
of  each  time  step  j#  as  calculated  from  equation  (2.19). 
Tabulated  data  of  the  curves  shown  in  figure  8  are  found 
in  Appendix  J.  Fireball  rise  was  not  considered  when  cal¬ 
culating  these  temperature  curves  because  the  case  of  sta¬ 
tionary  fireballs  is  less  complicated.  The  missile  being 
considered#  #41#  was  launched  at  the  same  time  that  silo  #1 
was  hit#  and  each  RV  landed  at  its  designated  ground  zero 
(a  1-cell  CEP) , 

The  cumulative  effect  of  thermal  radiation  from 
more  than  one  fireball  is  illustrated  by  plotting  the 
temperature  rise  for  each  burst  scenario.  Thus#  the  bottom 
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curve  represents  the  temperature  rise  due  to  one  burst, 
the  next  curve  represents  the  rise  due  to  two  bursts,  and 
so  on.  The  first  burst  occurred  at  t  ■  0,  and  the  times 
of  the  three  remaining  bursts  are  indicated  by  arrows. 

The  figure  shows  that  the  temperature  rise  for  a  series  of 
n  bursts  (0  <  n  £  4)  follows  the  temperature  rise  for  a 
series  of  n-1  bursts  until  tho  time  when  the  nth  burst 
occurs.  At  the  end  of  that  time  Btep,  the  temperature  is 
higher  than  the  temperature  from  n-1  bursts  because  of 
exposure  to  the  radiation  from  the  additional  burst.  Thus, 
the  curve  representing  the  temperature  rise  due  to  four 
bursts  is  not  smooth  as  it  increases  but  has  slight  irregu¬ 
larities  that  indicate  the  detonation  of  another  burst. 

Also,  the  four-burst  curve  decreases  mere  rapidly  from  its 
maximum  than  the  other  three  curves  because  the  contribu¬ 
tion  of  thermal  radiation  from  bursts  1  and  2,  measured  by 
k 

the  term  Aq  ,  is  small  or  zero  at  late  times.  In  addition, 
convective  cooling  is  more  effective  because  the  missile 
skin  i.s  so  hot. 

The  Effect  of  Fireball  Rise.  The  curves  in  figure  8 
were  calculated  for  the  case  of  fireballs  that  remained  on 
the  ground  at  ground  zero.  In  reality,  a  fireball  rises 
along  with  an  expanding  dust  cloud.  For  this  thesis,  the 
rising  fireball  was  assumed  to  be  located  at  the  top  of  the 
dust  cloud  and  centered  over  ground  zero.  Therefore,  the 
height  of  the  fireball  is  equal  to  the  height  of  the  dust 
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cloud  as  calculated  In  equation  (2.20).  When  fireball 
rise  is  considered,  the  slant  range  from  burst  to  missile 
is  calculated  from  the  fireball's  height  rather  than  from 
ground  zero . 

Table  IV  compares  the  temperature  rise  of  missile 
#41  for  both  stationary  and  rising  fireballs  in  a  four 
burst  scenario.  The  data  shows  that  temperatures  calcu¬ 
lated  with  fireball  rise  are  generally  higher.  The  differ¬ 
ence  in  the  temperature  of  interest,  the  maximum  tempera¬ 
ture  (at  j  ■  8),  is  '“S  ®K.  In  terms  of  probability  of 
damage,  £>d(704.5  °K)  »  .47  and  E>d(709.6  °K)  ■  52  (see  Appen¬ 
dix  G  for  an  explanation  of  how  to  calculate  P^(I)).  Thus, 
for  missile  #41,  the  case  of  fireball  rise  resulted  in  a 
slightly  higher  probability  of  damage.  However,  the  differ¬ 
ence  in  maximum  temperatures  would  not  be  readily  seen  if 
the  values  given  in  Table  IV  for  rising  fireballs  were 
plotted  in  figure  8. 

A  comparison  of  maximum  temperatures  for  other 
missiles  is  given  in  table  V.  The  data  shows  that  the  dif¬ 
ference  'between  maximum  temperatures  is  smaller  for  missiles 
launched  farther  away  from  the  bursts.  However,  since  the 
case  of  fireball  rise  is  more  realistic,  the  remaining 
results  presented  in  this  chapter  will  be  those  calculated 
when  considering  that  case. 


TABLE  IV 


TEMPERATURE  RISE  CONSIDERING  STATIONARY 
AND  RISING  FIREBALLS 


Missile  #41 

Launch  timet  0  eec 

Time  of  first  burst t  0  sec 

1-cell  CEP 

Time  Step 

Stationary  Fireballs 

Rising  Fireballs 

j 

t2  (°K) 

t2  (°K) 

1 

325.4 

325.4 

2 

277.5 

377.5 

3 

447.1 

446.8 

4 

505.7 

505.6 

5 

574.5 

575.0 

6 

638.7 

640.4 

7 

693.6 

697.1 

8 

704.5 

709.6 

9 

695.7 

701.9 

10 

677.8 

685.2 

11 

652.3 

660.3 

12 

622.0 

630.2 

13 

589.7 

597.6 

MiCttMUM 

TABLE  V 

TEMPERATURES  CONSIDERING  STATIONARY 

AND  RISING  FIREBALLS 

Missile 
Time  of 

launch  timet  0  sec 
first  burstt  0  sec 
1-cell  CEP 

Stationary  Fireballs 

Rising  Fireballs 

Missile 

# 

MaxT  (°K) 

MaxT  ( °K) 

21 

2561.6 

2615.1 

31 

1161.3 

1175.2 

41 

704.5 

709.6 

51 

511.5 

513.7 

61 

416.9 

418.0 

71 

366.0 

366.6 

81 

336.9 

337,2 
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The  Effect  of  Halving  the  Tima  Step*  As  mentioned 
in  Chapter  II,  the  finite  time  step  used  to  solve  the  dif¬ 
ferential  equation  for  Tj  was  chosen  to  be  tmax,  the  time 
of  the  second  thermal  maximum.  However,  it  is  possible 
that  a  smaller  time  step  could  give  more  accurate  results. 

To  test  this  possibility,  the  time  step  was  decreased  to 

t  mv/2.  Figure  9  shows  that  the  temperature  rise  oalou- 
in  ax 

lated  using  the  reduced  time  step  does  not  differ  much  from 
the  temperature  rise  calculated  using  tmax  as  the  time  step 
(see  Appendix  J  for  tabulated  data  of  figurr.  9).  More 
importantly,  the  maximum  temperature  reached  in  both  oases 
is  nearly  the  same  (~709  °K) .  Therefore,  in  order  to  mini¬ 
mize  the  number  of  necessary  calculations,  the  finite  time 
step  remained  as  tmax. 

Comparing  the  Cumulative  to  the 
Nonoumulatlve  Case 

In  the  cumulative  case,  a  series  of  four  bursts  is 
treated  as  one  event,  and  'one  ourve  is  plotted  to  lilies-  . 
trate  the  temperature  rise  of  the  missile  skin.  In  the 
noncumulativo  case,  a  series  of  four  bursts  is  treated  as 
four  independent  events  and  therefore  four  curves  are  needed 
to  illustrate  the  temperature  rise  due  to  each  of  the  bursts. 

The  difference  between  considering  four  bursts  as 
one  event  as  opposed  to  four  separate  events  is  shown  for 
missile  #41  in  figure  10  (see  Appendix  J  for  tabulated  data) . 
The  uppermost  curve  represents  the  cumulative  temperature 


44 


Missile  #41 


Fig.  9.  Comparing  Tiae  Steps  t  vs.  t  ,_/2  (1-cell  CEP) 


Missile  #41 


1 

rise  due  to  four  bursts,  while  the  four  bottom  curves  repre¬ 
sent  the  temperature  rise  due  to  eaoh  burst  separately.  In 
both  oases,  the  bursts  are  considered  to  have  a  1-cell  OEP. 

For  the  noncumulative  case,  the  highest  temperature 
reached  is  caused  by  burst  #1,  since  the  missile  at  the 
beginning  of  its  flight  is  closest  to  this  burst.  However, 
this  maximum  temperature  of  ~407  °K  is  much  lower  than  the 
maximum  temperature  of  *'710  °K  readied  in  the  cumulative 
case.  Furthermore,  the  probability  of  damage  for  missile 
#41  in  the  noncumulative  case  is  0,  while  for  the  cumulative 
case  the  probability  of  damage  is  **.32.  Assuming  that  the 
cumulative  case  is  correct,  calculating  the  temperature 
rise  by  considering  each  burst  to  be  on  independent  event 
severely  overestimates  the  missile's  chance  of  survival. 

The  next  section  further  emphasizes  this  difference  by  pre¬ 
senting  the  results  using  a  10-cell  CEP  for  each  burst. 

Results  Using  a  10-cell  CEP 

As  explained  in  Chapter  II,  a  10-cell  CEP  accounts 
for  any  possible  landing  position  that  an  RV  can  have.  Thus, 
the  probability  of  survival  calculated  using  a  10-cell  CEP 
is  a  more  realistic  number  than  that  calculated  uning  a 
1-cell  CEP  unless  the  value  of  CEP  is  very  small.  Tables  VI 
and  VIII  present  the  probabilities  of  survival  for  certain 
missiles  in  both  the  cumulative  and  noncumulative  cases. 

For  completeness,  values  are  given  for  both  rising  and  sta¬ 
tionary  fireballs. 
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TABLE  VI 


PROBABILITY  OP  SURVIVAL— 
NONCUMULATIVE  CASE 


Miaaila  launch  timat  0 
Tima  of  fir at  bur at t  0 
10-call  CEP 

1  DO 

s: 

Stationary  Piraballa 

Riaing  Piraballa 

Misaila  # 

pa 

r. 

25 

0 

0 

26 

.575 

.545 

27 

.588 

.555 

28 

.701 

.674 

29 

1 

1 

TABLE  VII 

PROBABILITY  OP  SURVIVAL- 
CUMULATIVE  CASE 


Miaaila  launch  timai  0 
Time  of  firat  burst*  0 
10-call  CEP 

aac 

aac 

Miaaila  # 

Stationary  Piraballa 

p. 

Riaing  Piraballa 

p. 

38 

0 

0 

39 

.0616 

.0481 

40 

.0645 

.0505 

41 

.523 

.480 

42 

.507 

.464 

43 

.538 

.496 

44 

.925 

.910 

45 

.927 

.913 

46 

1 

1 

For  this  thesis,  the  sure-safe  limit  is  defined  by 
a  probability  of  damage  of  .02,  and  the  sure-kill  limit  by 
a  probability  of  damage  of  .98.  Any  values  lower  than  .02 
are  rounded  to  0,  and  any  values  higher  than  .98  are  rounded 
to  1.  Thus,  the  values  of  0  and  1  that  appear  in  tables  VI 
and  VII  are  a  result  of  probabilities  of  damage  that  were 
1  and  0,  respectively.  Also,  any  missile  not  listed  on  the 
two  tables  has  a  probability  of  survival  equal  to  0  or  1, 
depending  on  the  missile's  silo  position. 

Another  way  of  presenting  the  valuos  given  in 
tables  VI  and  VII  is  by  using  figure  1  to  illustrate  the 
location  of  the  sure-safe  and  sure-kill  regions.  This  is 
done  in  figure  11,  with  no  distinction  made  for  stationary 
or  rising  fireballs  since  the  value  of  P§  will  not  change 
the  location  of  the  regions.  The  most  obvious  characteris¬ 
tic  of  figure  11  is  that  the  sure-kill  region  for  the  cumula¬ 
tive  case  extends  over  2500  m  further  downfield  than  the 
nonoumulative  case,  thus  killing  thirteen  more  missiles. 

Also,  the  region  between  the  sure-kill  and  sure-safe  limits 
is  larger  for  the  cumulative  case,  indicating  that  the 
probability  of  survival  baaed  on  distance  is  not  as  steep 
of  a  function,  as  shown  in  figure  12.  The  curves  drawn  in 
figure  12  Indicate  the  general  shape  of  P§  vs  ground  range 
but  are  not  actual  functions. 

The  results  of  the  10-aell  CEP  calculations  rein¬ 
force  those  of  the  1-cell  CEP  calculations.  Both  show  that 
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Fig-  12.  P  vs.  Dist 


treating  bursts  as  independent  events  does  not  adequately 
predict  the  temperature  rise  of  a  missile  akin  exposed  to 
thermal  radiation  from  a  series  of  bursts.  Failing  to  con¬ 
sider  the  cumulative  effect  will  underestimate  the  thermal 
threat,  which  could  be  an  expensive  oversight  in  terms  of 
nuclear  defense. 

Comparing  Cumulative  Thermal 
E f feats  to~¥IaiJ^Hects 

In  general,  blast  (overpressure)  effects  are  con¬ 
sidered  to  be  more  lethal  than  thermal  effects,  and  in  many 
oases  thermal  radiation  is  not  viewed  as  a  major  threat. 
Prior  to  this  thesis  research,  a  total  survivability  analy¬ 
sis  of  the  Dense  Pack  missile  system  was  performed  during 
NE  6.95,  Nuclear  Survivability  of  Systems.  The  analysis 
included  the  same  RV  walk  attack  as  this  thesis,  but 
required  a  random  launch  order  for  the  missile  field  and 
considered  all  bursts  as  independent  events.  The  blast 
aode  written  for  NE  6.95  included  a  "tail-chase"  of  the 
blast  wave  to  the  moving  target  and  used  sure-safe  and  sure- 
kill  overpressures  of  1.8  and  4.5  psi,  respectively.  In 
order  to  compare  cumulative  thermal  effects  with  (noncumulu- 
tive)  blast  effects,  this  blast  code  was  modified  so  that 
each  missile  was  launched  as  the  first  burst  occurred.  The 
code  was  then  run  to  determine  the  sure-safe  and  stare-kill 
regions  for  blast  effects,  and  the  results  are  given  in 
table  VIII. 
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TABLE  VIII 

RESULTS  OF  NONCUMULATIVE  BLAST  EFFECTS 


Missile  launch  timet  0  sec 
Time  of  first  burst:  0  sec 
10-cell  CEP 


Missile  # 


Probability  of  Survival 
0 

.0717 

.0956 

.590 

.626 

.750 

1 


Comparing  these  results  to  table  VI  shows  that 
noncumulative  thermal  and  blast  effects  have  nearly  the 
same  sure-safe  and  sure-region*.  However/  noncumulative 
blast  effects  are  overwhelmed  by  the  cumulative  thermal 
threat.  Although  this  may  not  be  the  case  if  cumulative 
blast  effects  are  considered/  the  results  do  indicate  that 
in  a  scenario  such  as  Dense  Pack,  thermal  radiation  is  more 
lethal  than  previously  expected  if  a  series  of  bursts  is 
treated  as  one  event. 


Observations  on  the  Cumulative  Case 

The  final  section  of  this  chapter  presents  two 
topics  that  were  not  discussed  in  the  objectives  of  the 
thesis  but  could  be  of  interest.  The  first  topic,  dust 
shielding,  is  of  interest  because  the  dust  entrained  in  the 
air  could  greatly  reduce  the  amount  of  thermal  radiation 


that  is  transmitted  to  the  missile.  The  second  topic  com¬ 
pares  the  difference  between  using  a  1-cell  CEP  on  a  10-cell 


CEP. 


The  Effects  of  Dust  Shielding.  A  contact  surface 


'3 
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burst  produces  a  large  amount  of  dust  and  debris  that  is 
carried  upward  with  the  rising  fireball.  This  cloud  will 
impede  the  transmission  of  thermal  radiation  to  a  target 
by  causing  additional  scattering  and  absorption  of  the  radi¬ 
ation.  The  dust  cloud  created  by  one  fireball  may  also 
decrease  or  even  block  the  thermal  radiation  emitted  by 
another  fireball. 

The  complex  behavior  of  dust  clouds  was  not 
included  in  this  study,  and  therefore  the  given  results 
would  tend  to  overestimate  the  temperature  rise  of  the 
missile  skin.  However,  an  attempt  was  made  to  determine  if 
dust  shielding  could  be  included  in  the  analysis  by  con¬ 
sidering  the  position  of  the  missile  and  the  fireball. 

The  situation  is  shown  in  figure  13. 

As  figure  13  indicates,  the  dust  cloud  is  mode  d 
as  a  cylindrical  shape  (McGahan,  1971:39-41).  Depending  on 
the  missile's  position,  the  missile  skin  will  be  subjected 
to  either  radiation  transmitted  through  the  atmosphere  or 
radiation  first  transmitted  through  a  thickness  of  dust 
cloud.  For  this  simplified  analysis,  the  transmission 
through  the  dust  cloud  was  considered  to  depend  on  the 
angle  <J>.  If  the  situation  was  such  that  the  missile  was 


Fig.  13.  Dust  Shielding 


above  the  top  of  the  dust  cloud  or  0  was  lass  than  20* , 
then  the  transmittanaa  t  would  be  the  value  calculated  for 
the  atmosphere  (see  Appendix  C) .  However,  if  $  were  greater 
than  20°,  then  t  -would  be  set  equal  to  0.  This  limitation 
only  applied  to  transmission  through  the  cloud  from  the 
burst  in  question,  not  to  transmission  through  a  cloud 
created  by  another  burst. 

Results  of  the  analysis  for  a  1-cell  CEP  are  given 


in  Appendix  L.  They  show  that  the  dust  clouds  have  no 
effect  on  transmittance  for  any  of  the  missiles  listed  in 
table  VII.  In  fact,  the  transmittance  is  only  affected 


for  missiles  whose  silos  are  within  ”2750  m  of  the  bursts 
(i.e.,  silo  #15  and  below),  but  the  decrease  in  temperature 
rise  is  not  enough  to  raise  the  missiles'  probability  of 
survival. 

Comparing  1-cell  CEP  and  10-cell  CEP  Results.  As 

explained  in  Chapter  II,  a  four  burst,  10-cell  CEP,  cumula- 

4 

tive  case  scenario  requires  10  calculations  of  maximum 
skin  temperature  to  determine  the  probability  of  survival 
of  a  single  missile.  Depending  on  the  computer  system  used, 
these  calculations  can  take  large  amounts  of  computer  time 
and  therefore  several  hours  may  be  needed  to  obtain  results 
for  all  of  the  missiles  of  interest.  In  contrast,  a  four 
burst  1-cell  CEP  scenario  requires  only  one  maximum  tempera¬ 
ture  to  calculate  the  probability  of  survival  for  a  single 
missile,  as  previously  shown  in  figure  8.  Thus,  it  would 
be  interesting  to  compare  1-cell  CEP  results  to  10-cell  CEP 
results  to  determine  if  a  1-cell  CEP  can  adequately  predict 
a  missile's  probability  of  survival.  These  results  are  pre¬ 
sented  in  table  IX  for  both  stationary  and  rising  fireballs. 
The  table  shows  that  a  1-cell  CEP  underestimates 

the  10-cell  CEP  P  for  missiles  #39  and  #40  and  overesti- 

0 

mater  p  for  the  remaining  missiles.  However,  the  differ- 
0 

ence  between  the  two  values  of  P  is  less  than  10%  for 

s 

missiles  #39  and  #40  and  less  than  1%  for  the  others.  In 
addition,  the  1-cell  CEP  results  do  not  change  the  cumula¬ 
tive  sure-kill  and  sure-safe  regions  shown  in  figure  11. 


TABLE  IX 


PROBABILITY  OF  SURVIVAL,  CUMULATIVE  CASE- 
1-CELL  CEP  VS  10-CELL  CEP 


Missile 
Time  of 

launch  timet 
first  burst! 

0  sec 

0  sec 

Stationary 

;  Fireballs 

Rising  Fireballs 

I 

■a 

.  e. 

Missile 

l-cellT 

10-cell 

l-cell 

10-cell 

# 

CEP 

CEP 

CEP 

CEP 

38 

0 

0 

0 

0 

39 

.0563 

.0616 

.0433 

.0481 

40 

.0592 

.0645 

.0456 

.0505 

41 

.527 

.523 

.483 

.480 

42 

.510 

.507 

.466 

.464 

43 

.543 

.538 

.499 

.496 

44 

.930 

.925 

.916 

.910 

45 

.933 

.927 

.919 

.913 

46 

1 

1 

1 

1 

Thus,  in  the  interest  of  reducing  computer  time,  a  1-cell 
CEP  scenario  can  be  used  to  calculate  the  probability  of 
survival.  The  resulting  decrease  in  computer  time  allows 
the  calculations  to  be  easily  performed  on  a  personal  com¬ 
puter  . 

The  above  conclusion  leads  to  the  consideration  of 
calculating  the  probability  of  survival  for  more  than  four 
bursts  using  a  1-cell  CEP.  According  to  Chapter  III  and 
Appendix  H,  the  maximum  number  of  bursts  that  will  affect 
the  temperature  rise  of  a  missile  skin  was  determined  to  be 
eight.  Table  X  presents  the  probability  of  survival  of  eight 
bursts  for  both  stationary  and  rising  fireballs. 
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TABLE  X 


PROBABILITY  OF  SURVIVAL,  CUMULATIVE  CASE  — 
EIGHT  BURSTS,  1-CELL  CEP 


Missile  launch  timet  0  sea 
Tima  of  first  bursts  0  sec 


1 

Missile  # 

Stationary  Fireballs 

p. 

Rising  Fireballs 

p. 

43 

0 

0 

44 

.0503 

.0243 

45 

.0516 

.0251 

46 

.452 

.341 

47 

.437 

.327 

48 

.460 

.349 

49 

.884 

.827 

50 

.886 

•  .830 

|  51 

1 

1 

The  results  are  also  shown  in  figure  14,  which  com¬ 
pares  the  sura-kill  and  sure-safe  regions  for  four  and  eight 
bursts.  The  sure-kill  region  for  eight  bursts  is  extended 
~ 1000  motors  down-range,  killing  an  additional  five  missiles. 
Thus,  the  consideration  of  eight  bursts  further  Illustrates 
the  fact  that  cumulative  thermal  effects  are  much  more 
lethal  than  noncumulative  thermal  and  blast  effects. 
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Fig.  14.  Sure-safe  and  Sure-kill  Regions, 
4  and  8  Bursts  (1-cell  CEP) 
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Chapter  V.  Conclusions  and  ReGonvmandationa 

Conclusion! 

Based  on  the  results  given  in  Chapter  IV,  the 
following  conclusions  are  drawn s 

1.  The  calculated  value  of  the  maximum  skin  tempera¬ 
ture  for  a  missile  exposed  to  a  series  of  bursts  depends  on 
the  way  in  which  the  bursts  are  treated.  Considering  a 
series  of  bursts  as  one  event  (cumulative)  gives  a  higher 
temperature  than  considering  each  burst  as  an  independent 
event  (noncumulative) .  A  higher  temperature  yields  a  lower 
probability  of  survival  for  the  missile. 

2.  Maximum  temperatures  calculated  when  consider¬ 
ing  fireball  rise  are  higher  than  those  calculated  for  sta¬ 
tionary  fireballs.  However,  the  difference  between  the  two 
cases  is  not  large  enough  to  distinguish  between  sure-safe 
and  sure-kill  regions  for  either  case.  Therefore,  fireball 
rise  as  treated  in  this  thesis  does  not  have  a  significant 
effect  on  the  results. 

3.  The  difference  in  maximum  temperatures  between 
the  cumulative  and  noncumulative  cases  extends  the  sure-kill 
region  of  the  cumulative  case  further  downfield,  making  the 
cumulative  case  more  lethal.  Specifically,  for  a  four 
burst  scenario  with  each  missile  launched  as  the  first 


I 

burst  detonates ,  the  cumulative  case  kills  thirteen  more 
missiles  than  the  nonoumulative  case. 

,  4.  Cumulative  thermal  ef facts  are  also  shown  to  be 

i  more  lethal  than  blast  effects,  which  contradicts  the  cur¬ 

rent  opinion  that  blast  is  the  primary  kill  mechanism. 

Based  on  this  conclusion  and  conclusion  #3,  the  failure  to 
consider  the  cumulative  case  will  underestimate  the  thermal 
threat  and  could  be  an  expensive  oversight  in  terms  of 
nuclear  defense. 

5.  The  results  given  are  considered  to  be  conserva¬ 
tive  because  the  analysis  did  not  account  for  any  thermal 
radiation  shielding  by  dust.  However,  dust  shielding  cannot 
be  adequately  modeled  using  simple  geometry  between  the 
missile  and  the  dust  cloud. 

6.  In  order  to  reduce  computer  time,  a  1-cell  CEP 
rather  than  a  10-cell  CEP  can  be  used  to  calculate  the  proba 
bility  of  survival.  If  this  is  done,  the  maximum  number  of 
bursts  considered  can  be  raised  to  eight.  This  results  in 
five  more  missiles  being  killed  and  further  emphasizes  the 
importance  of  considering  the  cumulative  case. 

Recommendations 

Baaed  on  the  assumptions  presented  in  Chapter  I  and 
the  observations  made  during  this  study,  the  following  recom 
mendations  are  made  for  further  study: 
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1.  Several  assumptions  made  during  the  study 
ignored  the  effect  of  dust  shielding  on  the  amount  of  radia¬ 
tion  transmitted  to  the  target.  The  expression  derived 

for  t  was  based  on  data  given  for  a  clear  atmosphere  and 
therefore  overestimated  the  transmittance.  The  curve  CT 
(figure  4)  used  to  determine  the  fraction  of  thermal  radia¬ 
tion  emitted  at  the  time  t  was  derived  for  an  air  burst , 
and  would  be  lower  for  a  surface  burst.  Finally ,  the 
assumption  that  the  fireball  is  located  at  the  top  of  the 
dust  cloud  and  the  simplified  geometric  analysis  of  dust 
shielding  did  not  consider  any  of  the  detailed  physics 
involved  in  dust  cloud  formation.  Since  the  amount  of  dust 
lifted  into  the  air  by  a  series  of  2  MT  bursts  could  be 
quite  large,  the  effect  of  the  dust  on  the  transmittance  of 
thermal  radiation  should  be  investigated. 

2.  The  fireball  was  assumed  to  be  an  isotropic 
point  source.  This  is  a  good  assumption  for  distances  that 
are  much  greater  than  the  fireball  diameter.  In  this  study, 
the  length  of  the  missile  field  was  20,300  m.  However,  the 
maximum  fireball  diameter  for  a  2  MT  contact  surface  burst 
is  expected  to  be  3700  m  (Glasstone  and  Dolan,  1977i71). 
Thus,  '.he  distances  between  burst  and  target  are  less  than 
five  times  the  fireball  diameter.  It  should  be  determined 
if  the  distances  are  great  enough  to  justify  the  assumption 
of  a  point  source.  If  not,  the  effects  of  a  finite  source 
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on  the  results  presented  in  Chapter  IV  should  also  be 
investigated, 

3.  This  study  did  not  consider  the  synergistic 
effects  of  heat  and  blast.  It  is  possible  that  the  increase 
in  temperature  would  sufficiently  weaken  the  missile's  skin, 
so  that  a  lower  overpressure  would  cause  the  same  amount  of 
damage  that  a  higher  overpressure  caused  when  considered 
alone.  Synergistic  ef feces  should  be  studied  to  determine 
if  thermal  radiation  poses  an  additional  threat  by  making  a 
system  more  vulnerable  to  blast  effects. 

4.  The  program  written  for  this  study  is  limited 
because  it  was  written  for  a  specific  missile  system  and  a 
four  burst  scenario.  A  more  general  program  would  be  easier 
to  adapt  to  other  situations  and/or  materials. 


Appendix  A.  An  Exp ra salon  for  the  Thermal 

Diffusion  Time 


In  Chapter  XX,  the  expression  for  the  time  required 
for  heat  to  be  conducted  through  a  slab  of  thickness  d  with 
thermal  diffusivity  K  was  given  as: 

tdiff  -  d2 /K  (A. 1) 


This  expression  can  be  derived  from  the  one-dimensional 
heat  transfer  equation  (Holman,  1972:102) 


9 


*T(X,t) 


i  3T(X,t) 

k  ”Td — 


(A.2) 


using  the  method  of  dimensional  analysis. 

As  the  name  implies,  dimensional  analysis  involves 
the  algebraic  manipulation  of  the  dimensions  of  physical 
quantities  to  provide  information  about  the  physical  pro¬ 
cesses  involved.  The  basis  of  dimensional  analysis  is  the 
simple  principle  of  dimensional  homogeneity,  which  states 
that  an  equation  is  complete  only  if  the  dimensions  of  each 
term  are  the  same,  or  homogeneous  (Parkhurat,  1964:16) . 
Equation  (A.2)  can  be  shown  to  be  complete  by  noting  that 
Infinitesimals  possess  the  dimensions  of  the  physical  ele¬ 
ments  they  represent;  (Parkhurot,  1964:21)  .  Therefore,  each 


1 
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expression  on  either  side  of  equation  (A. 2)  has  the  diman- 
2 

sions  of  6/l  ,  where  e  indicates  a  temperature  unit  and  L 
indicates  a  unit  length. 

An  important  application  of  the  principle  of 
homogeneity  is  the  derivation  of  an  unknown  quantity  based 
on  that  quantity's  functional  dependence  on  the  physical 
parameters  involved  in  the  problem  (Parkhurst,  1964:16). 
This  derivation  is  commonly  called  the  indioial  method,  and 
i.  illustrated  here  by  derivin9  an  axpreealon  tor 
from  the  following  terms  appearing  in  equation  (A. 2): 


, 

fundamental 

term 

description 

unite 

3T  ■  >AT 

temperature  differential 

0 

3x  -  >Ax 

distance  differential 

L 

K 

thermal  diffusivity 

L2/t 

3t  *  >At  «  > ‘td±f  f 

diffusion  tjme 

t 

The  first  step  in  the  indioial  method  is  to  write 
the  unknown  quantity  as  a  product  of  terms  to  unspecified 
powers  or  indices: 

[tdiff]  -  (ATa  Axb  K0]  (A. 3) 

The  dimensions  of  each  term  are  then  substituted  into  the 
above  equation  to  obtain: 

[t]  -  [ (6) a  (L)b  (L2/t)c] 
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or  [tl  »  [{0)a  (L)b+2c  (t)“C]  (A. 4) 

The  indices  of  each  fundamental  unit  on  the  left-  and  right- 
hand  side  are  then  equated  to  obtain  j  simultaneous  equa- 

tions  ,  one  for  each  unit: 

‘  *  ‘  .1  '  i»i.  ;  >  . 

0  *  a 

0  ■  b  +  2c 
1  *  -c 

i  ,  '  • 

I  ; 

Xn  general,  for  n  indices,  if  k  of  the  j  equations  are 

(  \  *  ,  ' ' 

independent,  then  any  k  of  the .indices  can be  solved  for 
ih  terms  of  n-k  other  indices.  For  this  case,  all  of  the 

1 '  i  •  »  ' .  :  '  1  • 

indices,  can  he  solved  fox  to  obtain t 

1  \  '■ 

a  **  0 

■,  »•  w  , ,  •» 

b  »  2 

c  »  -1  1  ' 

Therefore,  the  expression  far  t^^  can  be  found  by  sub¬ 
stituting  the  above  values  for  a,  b  and  c  into  equation 
(A. 3)  to  obtain i 

kdi££  '  4*?'  X'1 
or,  sinco  Ax  - 

tdiff  -  d2/K  (A. S) 
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Appendix  3 .  Missile  Characteristics 

The  missile's  velocity  V,  altitude  a,  down-range- 
distahce  (from  the  first  row  of  missiles)  drd,  and  flight 

i  .  ,  • 

path  angle  0  must  be  known  at  a  given  time  t  in  order  to 

■i-  '  ,  ■ 

calculate  the  slant  range  SR,  correction  factor  CP,  trans¬ 
mittance  t,  ambient  air  properties,  and  the  heat  transfer 
coefficient  h,  The  missile  characteristics  are  shown  as  a 
fun a t ion  of  time  in  figure  3.  In  this  figure,  the  down- 
range-distance  curve  refers  only  to  the  ground  distance  from 
the  missile  silo,  not  from  the  first  row  of  missiles. 

In  order  to  use  the  information  on  the  graph  in  a 
computer  program,  each  curve  was  reduced  to  a  data  file  of 
values  read  from  the  curve  at  every  second.  Values  were 
taken  from  t  *  0  to  t  ■  50  seconds  for  the  velocity,  alti¬ 
tude,  and  down-range-distanao  data  files,  and  from  t  *  0 
to  t:  »  70  sec  for  the  flight  path  angle  data  file.  The 
contents  of  these  files  are  presented  in  table  B-I.  At 
the  beginning  of  the  computer  program,  the  data  files  are 
read  into  the  arrays  VDATA(Sl),  2DATA(51)  ,  XDATA(51) ,  and 
ANGDAT  (71)  . 

The  missile's  velocity,  altitude,  and  down-range- 
distance  at  any  time  t  can  now  be  determined  in  one  of  two 


ways: 
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1.  if  t  ts  <_  50  seconds,  then  v,  z,  or  drd  can  be 
obtained  by  linear  interpolation,  where: 

t-^ctctj,;  t1  and  t2  are  whole  numbers  and 


t2~tl  "  1 

yl<y<y2;  ySv#  Z>  °r  drd 

and  y  ■  (t-t^) (yj-y^  +  yl  (B.l) 

2.  If  t  >  50  seconds,  then 

v  -  10 5t  -  1,350  (B.2) 

z  «  2 , 460 t  -  65,000  (B.3) 

drd  -  3,700t  -  127,000  (B.4) 


Since  the  curves  in  the  figure  are  in  terms  of  feet,  all 
parameter  values  must  be  converted  to  meters.  Also,  the 
values  of  drd  must  be  modified  to  reflect  the  missile's 
distance  from  x  -  0,  or  the  first  row  of  missiles.  This 
is  accomplished  by  adding  the  missile  silo's  x-coordinate 
to  drd: 

drd  -  drd*. 0348  4-  Sv(Mn)  (B.5) 

#* 

where  Sx(Mn)  -  Silo  x-coordinate  of  missile  Mn. 

The  missile's  flight-path-angle  0  is  found  in  a 
similar  manner: 

1.  If  t  <  70  seconds,  then 

8  -  (t-t^  (02-9;l)  +  ex  (B.  6) 
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2.  If  t  >  70  seconds ,  then: 

e  -  30.6  -  . 1 (t-70)  (B.7) 

Equation  (B.7)  is  a  linear  extension  of  the  curve  in 
figure  3  past  t  -  70  seconds.  The  value  of  6  is  then  con¬ 
verted  to  radians. 
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Appendix  C.  Calculating  the  Slant  Range ,  Correction 

Factor ,  and  Transmittance 

The  calculation  of  SR,  CF,  and  t  are  presented  in 
one  appendix  because  the  latter  two  terms  depend  on  the 
value  obtained  for  the  slant  range.  However,  the  procedures 
used  to  calculate  each  term  are  different  and  will  be  pre¬ 
sented  in  separate  sections  within  this  appendix. 

Calculating  the  Slant  Range 

The  slant  range  is  generally  defined  as  the  dis¬ 
tance  from  the  burst  to  the  target.  In  this  report,  the 
slant  range  is  the  distance  from  the  fireball  point  source 
to  the  location  of  the  missile  as  defined  in  figure  3. 

If  fireball  rise  is  being  considered,  the  expression  for 
SR  must  take  the  changing  position  of  the  point  Bource  into 
account.  The  slant  range  for  any  burst-missile  combination 
is  shown  in  figure  C-l.  The  missile  has  no  velocity  in  the 
y  direction  and  thus  travels  straight  north  along  the 
x-axis.  According  to  figure  C-li 

SR  -  [GR2  +  (alt  -  Hfb)2]0,5  (C.l) 

Hfb  is  calculated  using  the  VORDUM  equation  for  the  top  of 
the  dust  cloud  (McGahan  et  al.,  1971s40): 

Hfb  -  21 , G40 , 8w’ 177  [l-(l-t/t  )2]  (C.2) 
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Fig.  o-l.  Slant  Rang*  from  Bur at  to  Targat 

whara  w  -  yield  in  megatons,  and  tf  -  240a  (Bridgman,  1984), 
tha  tima  of  cloud  stabilization. 

For  an  RV  that  landa  at  ita  daaignatad  ground  zero 
(tha  targatad  ailo) ,  tha  axpraaaion  for  tha  ground  zaro 
range  ias 


GR  -  ((Mny  -  DGZy)2  +  (drd  -  DGZX)2) 


(C.  3) 


However,  if  tha  RV  haa  an  aiming  error  aaaooiatad  with  ita 
landing  poaition,  than  tha  axpraaaion  for  GR  muat  account 
for  tha  diaplacement  from  daaignatad  ground  aero.  Thia  dis- 
placamant  is  quantified  using  tha  10 -call  CEP  described  in 
Appendix  J .  The  situation  for  calculating  GR  for  a  call 
ia  shown  in  figure  C-2. 
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Fig.  02.  Ground  Range  from  Cell  i 


From  figure  C-2,  a  more  general  expression  for  ground  range 
is  given  by: 

-  .  5 

GR  -  ((Mny  -  yb)2  4-  (drd  -  x^2]*  (C.4) 

where  xb  -  DGZX  +  cos  8,^  and  yb  *  DGZy  +  sin  0^ 

For  the  middle  cell  (a  direct  hit),  equation  (C.4)  reduces 
to  equation  (C.3). 


I 
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Calculating  the  Correction  Factor 

In  Chapter  II ,  the  correction  factor  CF  was  defined 
as  sin  4/  where  t|>  is  the  angle  between  the  missile  skin 
surface  and  the  slant  range  vector.  The  correction  factor 
is  needed  to  calculate  the  amount  of  thermal  radiation  that 
is  incident  perpendicular  to  the  missile  skin's  surface. 

CF  can  be  calculated  by  recognizing  that  cos  i|>  is  defined 
by  the  dot  product  of  the  slant  range  and  velocity  vectors 
as  shown  in  figure  C-3> 


8  y  SR 


Missile 

skin 

surface 


6  ■  flight  path 
angle 


Fig.  C-3.  How  to  Calculate  CF 


By  defintion  of  the  dot  product i 


cos  “ 


SR„v„.  +  SR  v 


SR*  v 


(C.  5) 


'tXVUWL'UUUu' 


"  '"vVL.’W  UV  ' LS i 'lAa 


where 


v  *  v  cca  6 
#» 

v_  *  v  sin  9 
2 

SR  *  GR  «  drd  ~  x. 

XX  o 

SR.  *  alt  -  Hfb 
2 

The  product  of  y  components  does  not  appear  in  equatibn 
(C.  5)  because  V  *■  b»  It  should  also  be  noted  that  nega¬ 
tive  values  of  SR.  are  possible  if  alt  <  Hfb.  Although  a 
negative  value  does  not  mean  much  physiaally,  it  is  neces¬ 
sary  to  obtain  the  correct  answer  geometrically.  Once 
cos  V  is  known ,  CF  can  be  determined  from: 

CF  *  sin'}'  -  (1  -  cos2 'I')’5  (C.6) 

Calculating  the  Transmittance 

The  transmittance  t  is  defined  as  the  fraction  of 
direct  and  scattered  radiation  that  is  transmitted  to  the 
target  (Glassfcone  and  Dolan,  1977:316).  Because  t  is  a 
complex  function  of  atmospheric  conditions  and  distance,  a 

i  » 

detailed  analysis  of  the  correct  values  of  t  for  the  condi¬ 
tions  of  the  Dense  Pack  scenario  was  not  conducted. 

Instead,  t  is  assumed  to  behave  as  shown  in  figure  C-4 
(Glasstone  and  Dolan,  1977:318).  The  limitations  of  this 
figure  are  as  follows: 

1.  The  target  is  such  that  scattered  radiation  is 
received  from  all  directions. 
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2.  The  target  is  on  the  ground. 

3.  Visibility  is  12  miles  (21  km). 

Figure  C-4  is  used  to  express  x  as  a  function  of  distance 
for  a  surface  burst  by  reading  where  the  constant  lines  of 
x  intersect  the  bottom  of  the  graph  (zero  altitude) . 

Table  C-I  summarizes  the  values  of  x  taken  from  figure  C-4. 


TABLE  C-I 

VALUES  OF  X  AT  HEIGHT  OF  BURST  «  0 
(Glasstone  and  Dolan,  1977:318) 


X 

Distance  D  (kft) 

Distance  D  (m) 

0.01 

168 

51,206 

0.1 

85 

25,908 

0.2 

58 

17,678 

0.3 

50 

15,240 

0.4 

38 

11,582 

0.5 

31 

9,449 

0.6 

20 

6,096 

0.7 

15 

4,572 

0.8 

9 

2,743 

0.9 

5 

1,524 

Note:  All  distances  have  an  uncer'.ainty  of  -2  kft. 


When  plotted  as  x  versus  distance,  the  points  in 
table  C-I  appear  to  follow  an  exponential  curve  that  could 
be  fit  by  an  equation  of  the  form: 


exp (f (D) ) 


(C.  7) 


where  f(D)  represents  a  function  of  distance  D.  By  writing 
equation  (C.7)  as  Ln(x)  -  f(D),'  and  using  a  linear  regression 
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■  »  *  r  ’  i  *  w  i 

,fc, 


program  (Heilborn,  1981s 16) ,  the  following  expression  was 
obtained: 

£(0)  -  -.02455  -  6.439X10"5D  -  1. 407X10~9D2  +  1. 792X10“14D3 

(C.8) 

The  points  in  table  C-I  and  the  curve  produced  by  substi¬ 
tuting  equation  (C.8)  into  equation  (C,7)  are  shown  in 
figure  C-5.  The  coefficients  in  equation  (C.8)  were 
rounded  to  four  significant  digits,  and  the  standard  error 
estimate  of  f(D)  (the  standard  deviation  of  ln(x)  about  the 
fitted  curve)  is  0.05968.  This  error  was  considered  to  be 
well  within  the  error  inherent  in  reading  from  figure  C-4. 

Although  the  expression  for  x  was  derived  from 
data  for  a  surface  burst  and  ground  target,  the  distance  D 
in  aquation  (C.8)  is  taken  to  be  the  slant  range.  This  was 
done  because  slant  range  represents  the  distance  in  the 
atmosphere  that  the  thermal  radiation  must  travel  to  reach 
the  missile.  The  expression  for  x  also  does  not  consider 
the  changing  atmospheric  conditions  as  the  fireball  rises 
or  the  effects  of  dust  shielding  the  fireball.  However, 
since  x  is  inherently  uncertain  because  of  its  complex 
dependence  on  many  variables,  the  expression  derived  for 
use  in  this  analysis  is  considered  to  be  adequate  enough  to 
show  the  change  in  transmittance  aB  the  missile  moves  away 
from  the  burst. 
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Appendix  D.  Ambient  Air  Properties 


As  mentioned  in  Chapter  II,  values  for  the  ambient 
air  temperature  T  ,  density  p  .  dynamic  viscosity  u  ,  and 

a  ft  a 

conductivity  k  at  the  missile's  altitude  are  needed  to 

cl 

calculate  the  heat  transfer  coefficient  and  skin  temperature. 
These  values  are  obtained  using  equations  given  in  the 
reference  from  the  National  Oceanic  and  Atmospheric  Admin¬ 
istration  (NOAA,  1976s 6-20). 

The  general  equations  for  air  temperature  and 
pressure  ares 

IP  Lki»0i  T(z)  -  Tk  +  Lk(z  -  zk) 

T.  ( . 034164/L.  ) 

p<d  -  Pk  (d,1) 

IF  Lk-  0:  T(z)  -  Tk 

- . 034164  ( z-z.) 

P(z)  -  Pkexp  [ - <5 - —  ]  (D .  2) 

k 

where  z  is  the  missile's  altitude  in  meters,  T(z)  is  the 

2 

temperature  in  °K,  and  P(z)  is  the  pressure  in  N/m  .  The 
values  of  Lk,  Tk,  and  Pk  depend  on  the  region  of  the 
atmosphere  and  are  given  in  table  D-l. 
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TABLE  D-I 


VALUES  USED  TO  CALCULATE  AMBIENT  AIR 
TEMPERATURE  AND  PRESSURE 


Region 

z  k  (m) 

Lk 

Tr  (°K) 

Pk  (N/m2) 

0  <  z  <  11,000 

0 

-0.006545 

288.15 

1.013X105 

11,000  <  z  <  20,000 

11,000 

0.0 

216.65 

2.269X104 

20,000  <  z  <  32,000 

20,000 

0.0010 

216.65 

5.528X103 

32,000  <  z  <  47,000 

32,000 

0.0028 

228.65 

8.888X102 

once  T ( 2)  and  P(2)  are  known,  the  remaining  terms 
are  calculated  as  follows: 


pa  -  .003484  !|fj-  kg/m3  (D.3) 

»  1.458X10"6[T(z)]1,S  kg/m-s  (D.4) 

ya  T  ( z ) +1 10 .4 


k  «  2.64638X10~3[T(z) l1,5 
a  T(z)+245.4(10 ^ -i2/T(z) ) j 


J/m-s~°K  (D.5) 
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Appendix  E.  Local  Leat  Transfer  Coefficient 

The  local  convective  heat  transfer  coefficient#  h# 
was  introduced  in  equation  (2.14)#  which  describes  the  rate 
of  energy  loss  due  to  convection  (Holman#  1976:12): 

^convection  "  -Tair  (E.l) 

2 

where  h  has  units  of  J/m  -s-°K 

The  expressions  for  h  used  in  this  report  were 
those  derived  for  viscous  flow  over  a  flat  plate.  The 
derivation  involved  a  detailed  mathematical  analysis  of 
the  thermal  boundary  layer#  the  region  where  temperature 
gradients  are  present  in  the  flow  along  the  plate  (Holman# 
1976:154-171).  The  results  depend  on  the  type  of  flow# 
which  is  described  by  the  Reynolds  number.  For  viscous 
flow#  the  Reynolds  number  is  (Holman#  1976:149) 

Re  -  (E.2) 

y 


where  v  »  fluid  free-stream  velocity  (m/s) 
p  -  fluid  density  (kg/m^) 

]x  *  fluid  dynamic  viscosity  (kg/m-s) 
x  *  distance  from  leading  edge  of  plate  (m) 

For  air,  p  =*  p.  and  y  «  u  (see  Appendix  D)  .  For 

cl  cl 

this  report,  the  velocity  was  taken  to  be  the  speed  of  the 
i 
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missile,  and  the  value  for  x  “  xm  was  chosen  to  be  5.5  m, 

as  explained  at  the  end  of  this  appendix.  For  flat  plates, 

the  transition  between  laminar  and  turbulent  flow  is  con- 

5 

sidered  to  occur  at  Re  »  5X10  (Holman,  1976:149). 

Another  term  used  to  calculate  the  heat  transfer 
coefficient  is  the  Prandtl  number,  which  expresses  the  rela¬ 
tive  magnitudes  of  momentum  and  heat  diffusion  in  the 
fluid.  The  expression  for  the  Prandtl  number  is  (Holman, 
1976:164) 

cnw 

Pr  -  -g-  (E.3) 


where  Cp  -  specific  heat  capacity  of  the  fluid  (J/kg-K) 
K  ■  thermal  conductivity  of  the  fluid  (J/s-m-K) 

For  air,  c  ■  1005  J/kg-K  (Holman,  1976:503),  y  ■  y  ,  and 
P  a 

k  ■  k^  (see  Appendix  D) . 

.From  the  analysis  of  the  thermal  boundary  layer, 
the  following  expression  is  derived  (Holman,  1976:170): 


Nu 


so  that 


h  -  ^ 


(E.4) 


(E.5) 


where  Nu  is  the  Nusselt  number.  For  a  plate  heated  over 
its  entire  length,  the  equations  for  Nu  are 
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(E.6) 


if  Re  <  5X10%  Nu  -  . 322Pr1/3Re1/2 
(Holman,  1976:171) 

if  Re  >  5X105:  Nu  -  . 0296Pr1/3Re4/5 

(Holman,  1976:180)  (2.7) 


Thus,  h  can  be  calculated  from  equation  (E.5)  if 
the  missile  velocity  and  several  air  properties  are  known. 
As  an  example,  at  t  *  5  seconds  into  the  missile's  flight, 
the  missile  velocity  v  ■  140  ft/sec  -  42.8  m/s  and 
z  -  600  ft  *  182.9  m.  The  following  values  of  ambient  air 
properties  at  z  •  182.9  m  can  be  calculated  from  the  equa¬ 
tions  in  Appendix  D: 

Pa  -  1.203  kg/m3 
u  «  1.784X10’5  kg/m-s 

ft 

k,  -  .02523  J/m-s°K) 

A 

Using  these  values  in  equations  (E.2)  and  (E.3)  results  in 


and 


Re  .  (42.8)  (5.5)  (1.203) 
(1.784X10"5) 


1.5 9X10 7 


or  (1005)  ( 1 . 784X10**5) 
Fr  (762523) 


.711 


5 

Since  Re  >  5X10  ,  the  flow  is  turbulent  and  the  Nusselt 
number  is  given  by  equation  (E.7): 
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NU  ■  .0296(.711)1/3(1.59X107)4/5 
-  1.52X104 

Substituting  this  value  of  Nu  into  equation  (E.2V  produces 

i 

a  value  for  the  heat  transfer  coefficient  h: 

h  -  (• 0*523^1. 5m0<|  .  69_7  J/m2.s.„K 

Choosing  a  Value  for  x^ 

Figure  E-l  shows  the  variation  of  h  with  Reynolds 

number  for  several  values  of  x  .  where  is  the  distance 

m  m 

from  the  leading  edge  of  the  missile.  The  data  shows  that 

the  choice  of  x^  *  5.5  m  represents  an  average  between  a 

maximum  amount  of  heat  transfer  occurring  for  x  »  1.0  m 

m 

and  a  minimum  amount  of  heat  transfer  occurring  at  xm 
■  20  m. 


« 
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Variation  of  h  with  Reynolds  Number  for  Various  Values  of 


Appendix  F.  Calculating  the  Fraction  of  Thermal 

Energy  Emitted,  ACT  j 

The  term  CT  refers  to  the  curve  shown  in  figure  4 
in  Chapter  II.  This  curve  represents  the  fraction  of 
thermal  energy  emitted  up  to  any  time  t'  *  where  t 

IUClX 

0  44 

is  the  time  after  detonation,  t  ■  0.417V  (Glasstone 

max 

and  Dolan,  1977:310) ,  and  Y  is  the  yield  in  kilotons.  At 

any  time  step  j  in  the  computer  calculation  of  the  missile 

) 

skin  temperature,  the  total  thermal  fluenae  that  is  incident 
on  the  surface  during  that  time  step  is 

l 

AQj  -  ACTj  •  Finflident  (*.l) 

where  F^ro^ent  is  ^i,ven  equation  (2.13).  The  curve 
in  figure  4  must  be  digitized  in  order  to  calculate  ACT^ 
during  the  run  of  the  computer  program,  Sections  of  the 
curve  wero  fit  with  second-degree  aquations  or  less  using 
points  from  the  curve  and  a  linear  regression  program 
(Heilborn,  1981:16).  Table  F-I  presents  the  points, 
equations,  and  correlation  coefficients  for  the  equations. 
Values  of  CT  for  t'  <  1.5  ware  read  directly  from  figure  4, 
while  the  points  for  t"  >_  1.5  were  obtained  with  the  help 
of  a  digitizer  (owned  by  the  Plasma  Physics  Group,  Aero 
Propulsion  Laboratory,  Bldg  450,  Wright-Patterson  AFB,  Ohio) 


Values  Used  Equation  (froae  Correlation 

Region  t "  CT  Linear  Regression)*  Coefficient** 


co 

co 
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TABLE  F-2 


DATA  FROM  DIGITIZER 


1.286297 

0.295778 

1.S28967 

0.367067 

1.771637 

0.426614 

2.014308 

0.470785 

2.292048 

0.507967 

2.534717 

0.532569 

2.815257 

0.560525 

3.057937 

0.579815 

3.300607 

0.598266 

3.543277 

0.615599 

3.785947 

0.632653 

4.066487 

0.649986 

4.309158 

0.659770 

4.620568 

0.678222 

4.863238 

0.687727 

5.105907 

0.698071 

5.441158 

0.709253 

5.822697 

0.721274 

6.184597 

0.733296 

6.427267 

0.736650 

6.669937 

0.744758 

6.912607 

0.750908 

7.155278 

0.756779 

7.397948 

0.762091 

7.640628 

0.769080 

7.883297 

0.770478 

8.125967 

0.776628 

8.368637 

0.780262 

8.611307 

0.783897 

8.853978 

0.789488 

9.096648 

0.791445 

+t'  -  *  -.0169327  ) 

>  u'  and  CT  normalized  to  0,0. 
CT  *  CT  *  +  .167738  5 
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f- 

Given  the  equations  in  table  F-I,  the  value  of  CT 
at  any  t'  can  be  calculated  to  obtain  ACTj.  As  an  example, 
for  the  second  time  step  ( j-2) ,  the  values  of  t'  at  the 
beginning  and  end  of  that  time  step  for  burst  #1  are  1  and 
2  respectively  (see  figure  7  in  Chapter  III) .  Using  these 
values  of  t' ,  ACT2  is  then: 

ACT2  -  CT  (t '*■2)  -  CT(t^-l) 

-  0.4680  -  0.20 

-  0.268 
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Appendix  G.  The  Probability  of  Damage 

Function  and  CEP 

The  Probability  Function 

This  appendix  contains  methods  for  calculating  the 
probability  of  damage  and  quantifying  the  aiming  error. 

The  information  was  taken  exclusively  from  class  notes  given 
by  Dr.  C.  Bridgman  (Bridgman,  1984) . 

For  this  report,  probabilities  of  damage  are  calcu¬ 
lated  using  intensity  I  as  the  independent  variable  and  the 
cumulative  log  normal  function  as  the  distribution  function. 
Each  intensity  has  an  associated  probability  of  damage 
according  tot 


pdm 


1 

/2tt0'I 


exp  [- 


1  ilnl-a'* 

I(  a' 


2 

)  ]  dl 


(G.l) 


where  a'  and  3"  are  parameters  that  depend  on  the  sure- 
safe  and  sure-kill  intensities  of  the  material.  This 
integral  can  be  solved  analytically  using  the  following  sub¬ 
stitution: 


u 


In  1  -a" 
3' 


(G.  2) 


Substituting  equation  (G.2)  into  the  integral  of  equa¬ 
tion  (G.l)  produces 


1_ 

2rr 


du 


(G.3) 


I 

N 

I 


P.(u)  mfU  — - —  e"u2/2 

Q  /  27T 


whoie  integrand  i»  recognized  as.  the  normal  distribution 
function  about  u  *  0.  A  graph  of  the  Above  integral  is 
shown  in  figure  G-l. 


1 

99 

1 

1 

K 

1 

I 


H 

I 


Equation  (G.3)  can  be  solved  in  any  number  of  ways. 
A  good  method  for  computer  implementation  is  to  use  the 
following  equations  (Abramowitz  and  Stegun,  1965t932)t 


P(u1-  1  - 


[2(l+.196854u'+.  11594 (u-)  2+.  000344  (u')  3+. 019527  (u') 4]  * 


(G.4) 
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If  u 


>  o,  Pd(I)  -  P(u>) 

If  u  <0,  Pd(I)  -  1  -  P(u') 

Thus,  once  an  intensity  I  is  found  and  and  0"  are  known, 
Pd(I)  can  be  calculated  using  equations  (G.2)  and  (G.4) 
or  equation  (G.2)  and  the  curve  in  figure  G-l. 

As  mentioned  before,  the  parameters  a'  and  are 

calculated  from  sure-safe  and  sure-kill  intensities  I  .  and 

88 

I&kr  For  analysis,  all  probabilities  of  damage  below 

0.02  are  considered  to  be  sure-safe  (Pd(I)»0),  and  all 
probabilities  above  0.98  are  considered  to  be  sure-kill 
(Pd ( I ) -1.0) .  Thus,  a'  and  0“  must  be  defined  so  that  if 
Pd(I§k)  ■  0.98  and  Pd<IBg)  ■  0.02  then  the  following  are 
trues 

f  8*  1  i  ini-a'  ^ 

.98-1  — - -  exp  [-  -4-  ]  dl  (G.  5) 

Jo  /T T  Q'  i  2  6 


.02  -f  88  - r -  exp  [--i-  (—A"  %')2]  dl  (G.  6) 

J0  /  2tr  0 '  I  2  p 


The  solution  to  the  above  equations  yields 


In  I  -  oc' 

3  S 


-2.054 


(G.7) 


In  I 


sk 

0' 


a 


2.054 


(G.  8) 
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which  can  in  turn  be  solved  to  give 


2  ^ss^sk* 

(G.9) 

TTOT  ln  (l  ] 

ss 

(G. 10) 

The  expressions  for  ol'  and  S'*  given  by  equations  (G.9) 
and  (G. 10)  are  used  to  calculate  u  in  equation  (G.2) 

Aiming  Error  or  CEP 

The  development  of  the  function  describing  the 
aiming  error  of  an  RV  around  designated  ground  zero  (DGZ) 
assumes  the  following! 

1.  Errors  are  in  the  x  and  y  directions  only.  The 
height  of  burst  (z«0)  is  accurate. 

2.  The  errors  in  x  and  y  are  entirely  random. 

3.  o  m  a  . 

x  y 

If  f(x)  *  f(y)  are  functions  describing  the  errors  in  x  and 
y  respectively,  then 

1  i  x  ^ 

f  (X)  -  — —  exp  [-i  (f-)  ]  (G.  11) 

✓2*0  2  x 

X 

1  1  2 

f  (y)  -  «XP  [-4  <?->  1  (G.  12) 

✓27c  *  y 

y 

These  two  distribution  functions  are  shown  in  figure  G-2 
as  two  gaussians  centered  at  the  aim  point  (DGZ)  and 


97 


extending  to  define  a  circular  area  A.  The  probability  of 
hitting  area  A  ia 

2  2 

f(x,y)  -  f(x)f(y)  ■ — exp[- -j  ]  (G.13) 

2 1T0  0* 


which  is  a  circular  area  density  or  planar  normal  function. 


Fig.  G-2.  Planar  Normal  Function 


Equation  (G.13)  can  be  written  for  a  differential  area  dA 
as 

1  i  r  2 

f(r)dA  -  — exp[-'j  (-)  ]  2wrdr  (G.14) 

2no 

2  2 

where  r  *  x  +  y  ,  the  radius  of  the  circle  defined  by 
f(x,y).  The  probability  of  falling  anywhere  inside  of 
radius  r  is  then 

F(r)  -fr  f (r)dA  -  f*  — ^  exp[-i(|)2  2wrdr  (G.15) 
JO  Jo  2wo2  2  0 


F(r)  is  a  cumulative  circular  normal  function.  By  sub¬ 
stituting  r  -  so  and  recognizing  that  £(s)  -  £(r)(dr/ds), 
the  above  equation  can  be  written  as 


r»a  i  2 

F (■)  -J  exp[--j*  J  sds 


(G. 16) 


F(r)  can  be  evaluated  for  different  values  of  s  as  follows i 


0.393 


0.865 

0.989 


The  term  circular  error  probable,  or  CEP,  is 
defined  as  the  radius  r  inside  of  which  501  of  the  tar¬ 
geted  RV's  will  hit.  According  to  the  above  values  of 
F,  this  radius  is  between  0  and  20  .  Specifically,  using 
equation  (G.15) 


f CEP  ,  ,2 

,50  ml  - 2  exp|;-^(~)  ]  2iTrdr 


(G. 17) 


By  substituting  w  ■  0.5(r/o)  ,  the  above  integral  can  be 


solved  analytically  to  obtain 


CEP  -  /2  In  2  ~  1.180 


(G. 18) 


This  result  will  be  used  to  define  areas  of  equal  proba¬ 
bility  around  the  designated  ground  zero. 


.V-.  „L.  . 
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Combining  Probability  of  Damage 


Lng  Error 
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The  overall  probability  of  damage  to  a  missile 
depends  on  the  following: 

a.  the  probability  of  the  RV  hitting  any  dA  sur¬ 
rounding  the  DGZ 

b.  the  probability  of  damage  given  a  hit  on  any  dA 

The  condition  in  (a)  is  given  by  equation  (G . 15 ) .  Rewriting 
equation  (G.15)  in  polar  coordinates ,  where  r  -  p  and 
dA  «  dpde ,  yields: 


f  (p)  dA  ■  a  expf-  x  (§)  ]  PdPdQ 
2’n‘C*  *  ° 


(G. 19) 


The  condition  in  (b)  is  given  by  Pd(I)  (equation  (G.l)). 
Thus,  the  probability  of  damage  over  all  dA  is 


n2TT 

(IP)  Pd(I)dR 

■f  expl-^f,2]/"1  —A -  )  dl  jpdpd0 

Jo  Jo  1 2^0  2  Jo  2  p  ‘ 


(G.  20) 


This  integral  can  be  simplified  to  a  sum  over  N  areas, 
designated  as  A  -  . ,  as  follows: 

©  ci  x 


?d  "  ^  f<»i  pdi  ^^area  i 


(G .  21) 


100 


I 


The  areas  are  now  chosen  so  that 


f(p)iA 


area  i 


(G. 22) 


That  is,  the  space  around  designated  ground  zero  is  divided 
into  areas,  or  cells,  of  equal  probability  of  getting  hit. 
Thus,  the  probability  of  damage  for  a  missile  becomes 


pd  “  §  2  pd  (I) 

i®l  i 


(G.  23) 


The  probability  of  damage  from  a  hit  in  cell  i,  P.  (I) , 

di 

is  calculated  by  finding  the  intensity  I  at  the  slant  range 

determined  from  the  centroid  of  cell  i  (see  Appendix  C) , 

and  then  using  equations  (G.2)  and  (G.4).  The  P.  (7.)  '  b 

aL 

found  for  each  cell  are  summed  and  divided  by  the  total 
number  of  cells,  N,  to  obtain  P^.  Knowing  P^,  the  proba¬ 
bility  of  survival  P#  ■  1  -  P^. 

Calculating  the  Cell  Centroids 

A  10-cell  space  around  DGZ  is  shown  in  figure  G-3. 
For  convenience,  each  cell  is  numbered.  The  dot  inside 
each  cell  is  the  cell  centroid,  whose  position  must  be 
known  to  calculate  the  slant  range.  Each  centroid  is 
defined  by  <P^>  and  <0^>,  where  <P^>  is  the  distance  to  the 
centroid  of  cell  i  from  DGZ  (the  center  of  cell  #1) ,  and 
<0£>  is  the  angle  between  <p^>  and  the  x-axis.  <p^> 


Fig.  G-3 .  A  10-call  Diacratizad  Spaca  Around  DGZ 

and  <6^>  art  darivad  auch  that  thay  daflna  cell*  of  aqual 
probability. 

A  nuxnbar  of  atapa  ara  naooaaary  to  obtain  an 
axpraaalon  for  <p^>  vialng  th«  following  variablaai 

p ^ *  outar  radiua  of  call  i 
nri  number  of  calls  in  ring  r 
Nr»  number  of  calla  inaida  ring  r 
:  total  number  of  calls 
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A  ring  of  cells  is  a  group  of  cells  with  equal  radii,  such 
as  cells  #2,  3,  4,  and  5.  Therefore,  for  any  cell  in  ring 
ri'  pi  of  that  is  synonymous  to  pr  for  the  ring.  For 

a  10 -cell  configuration,  the  following  table  can  be  con- 

r  i 

l  i 

4  5 

5  10 

Since  equal  probability  cells  are  desired,  the 
total  probability  of  hittinq  inside  ring  r  is  Nr/N^,  and 
in  turn  this  is  equal  to 


structedi 

■  Ring 
1 
2 

1  3 


(G. 24) 


for  all  6.  The  above  integral  can  be  solved  analytically 
using  the  substitution  w  -  p/a  to  obtain 


i  Pr  2 

U*  *  1  -  )  1  (G.  25) 


Solving  equation  (G.25)  for  pr  results  in 


P 


r 


o[-2  In  (1 


H 


(G. 26) 


Equation  (G.26)  can  be  generalized  by  expressing  a  in  terms 
of  CEP  using  equation  (G.lJ)t 
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Msammim 


N 


-JL 

CEP 


- 


~ lnU~  y 

In  2 


The  following  aquation  is  usad  to  find  <P  > i 


/'■ 


pf(p)  2irpdp 


^3?“1 
<P  >  -  L 

r  r  P, 


r r 

p  , f 

^r-l 


(p)  2irpdp 


(G. 27) 


(G. 28) 


Substituting  tha  axprassion  for  f(p),  and  using  equation 
(G.  24)  yields 

NT  f  1  P  2  D2 

<pr>  ,B  n  0  -»  1  %  dP  (G. 29) 

r  Pr-1  0 

Equation  (G.29)  is  solved  by  making  tha  substitution 
w  «  p/a  and  than  integrating  by  parts.  In  terms  of  CEP, 
tha  result  is 

«  • 

<Pr  >  Nm  (  P-.  Pr„l  2  p„  Pj.  2 

cep  "  tT  J  “ceF  ®Xp ln  2  ("ceF)  1  “  CEP  exP ln  2  {CEP>  1 

it  •  ®  •  S  P  -  .  5  P  j.  . 

+  <3~2>  [CNF  ( (2  ln  2)  ^)  -CNF((21n2)  -§£)] 

(G. 30) 

where  CNF(x)  »  cumulative  log  normal  function  for  an  argu¬ 
ment  x,  as  evaluated  in  equation  (G.4) .  For  the  10-cell 


104 


system,  the  following  values  are  obtained  from  equations 
(G. 27)  -  (G. 30) t 


TABLE  G-I 


VALUES  OP  <Pi> 


Ring  r 

Pi/CEP 

<Pi/CEP> 

1 

0.3899 

0.0 

2 

1.0 

0.7109 

3 

00 

1.509 

Thus,  <Pr>  can  be  found  by  multiplying 

<Pr/CEP>  by  the 

given  value  of  CEP. 

Table  G-II  shows  values  of  6^ 

found  by  inspection 

of  figure  G-3. 

TABLE  G-II 

VALUES  OP  <6^> 

Cell  #i 

ei 

1 

0 

2 

0.25" 

3 

0 . 751T 

4 

1.25" 

5 

1.75" 

6 

0 

7 

0.4" 

8 

0 . 8" 

9 

1.2" 

10 

1, 6" 

1 


Appendix  H>  Date  for  Choosing  the  Upper 

Limit  of  Bursts 


As  kentioned  in  Chapter  .III,  the  upper  limit  to 
the  number  of  bursts  that  will  cause's  missile's  akin 
temperature  to  rise  was  determined  to  be  eight.  This  limit 
was  estimated  by  calculating  the  maximum  temperature  reached 
for  a  series  of  bursts  and  assuming  that  each  RV  lands  at 
its  desighatad  ground  ?.ero  (a  1-cell  CEP) «  The  following 
table  summarises  the  data  ustfed  ho  arrive  at  the  upper  limit 
by  including  the  results  for  4,  7,  a  and  9  bursts.  It 
should  be  noted  that  missiles  below  #13  required  fewer 
than  eight  bursts  to  reach  a  maximum  temperature,  but 
these  missiles  have  no  probability  of  survival  for  even 
four  bursts. 


■  ^.yuv>vv/>^wuuuu:v.^\jui/wA:wu^  n««m  ,vu«.ymuh  uwiaaaafli  uwrauuuuafte 


TABLE  H-l 


i 

i 


DATA  USED  TO  DETERMINE  THE  UPPER  LIMIT  OF  BURSTS 

FOR  A  1-CELL  CEP 


Maximum  Temperature  Reached  ( °K)  * 


Missile 

# 

4  Bursts 
j-8 

7  Bursts 
j-12 

8  Bursts 
j-13 

9  Bursts 
j-13 

21 

2,615 

3,602 

3,611 

3,611 

23 

2,576 

3,543 

3,564 

3,564 

31 

1,175 

1,550 

1,565 

1,565 

33 

1,169 

1,541 

1,558 

1,558 

37 

891 

1,145 

1,157 

1,157 

41 

710 

886 

895 

895 

43 

708 

883 

893 

893 

47 

594 

721 

729 

729 

51 

514 

607 

612 

612 

53 

513 

606 

611 

611 

57 

458 

528 

532 

532 

61 

418 

471 

474 

474 

63 

418 

470 

473 

473 

67 

389 

429 

431 

431 

71 

367 

398 

399 

399 

73 

366 

397 

399 

399 

77 

350 

374 

375 

375 

81 

337 

356 

357 

357 

83 

337 

356 

357 

357 

93 

320 

331 

331 

331 

*Temperatures  are  rounded  to  the  nearest  degree; 
however ,  all  maximum  temperatures  for  8  and  9  burst*  were 
equal,  before  rounding. 

j  “  timeatap  when  maximum  temperature  was  reached. 

Fireball  rise  considered  for  all  results. 

Missiles  were  launched  when  silo  #1  was  hit  (tb, 

»  0  sec) .  1 
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Appendix  I.  Data  for  Choosing  the  Maximum 

Number  of  Time  Steps 

a 

As  mentioned  in  Chapter  III,  the  maximum  number  of 
time  steps  needed  to  determine  the  maximum  skin  temperature 
was  chosen  to  be  11.  This  choice  was  based  on  the  results 
from  a  1-cell  CEP,  4  burst  scenario  for  several  missiles. 
Table  1-1  summarizes  the  results  for  eight  missiles  along 
the  length  of  the  field. 


were  launched  when  silo  fl  was  hit.  (tb1  =  0  sec) 


Appendix  J.  Data  for  Figures  in  Chapter  IV 


This  appendix  contains  tables  of  data  that  were 
used  to  plot  figures  8,  9,  and  10. 


TABLE  J-I 


DATA  USED  TO  PLOT  FIGURE  8,  TEMPERATURE  RISE  FOR 
THE  CUMULATIVE  CASE  (1-CELL  CEP) 


1 

Burst 

2 

Bursts 

3 

Bursts 

4 

Bursts 

j 

t2(°k) 

j 

T2(°K) 

j 

t2(°k) 

j 

t2(°k) 

1 

325.40 

1 

325.40 

1 

325.40 

1 

325.40 

2 

374.45 

2 

377.49 

2 

377.49 

2 

377.49 

3 

392.67 

3 

447.08 

3 

447.08 

3 

447.08 

4 

400.78 

4 

491.91 

4 

505.69 

4 

505.69 

5 

405.62 

5 

508.35 

5 

574.46 

5 

574.46 

6 

406.99 

6 

514.60 

6 

604.06 

6 

638.73 

7 

404.98 

7 

513.99 

7 

608.34 

7 

693.58 

8 

400.28 

8 

507.85 

8 

604.07 

8 

704.54 

9 

393.33 

9 

496.97 

9 

592.34 

9 

695.66 

10 

385.51 

10 

483.00 

10 

575.02 

10 

677.84 

11 

466.01 

11 

552.79 

11 

652.34 

12 

527.99 

12 

622.02 

13 

502.76 

13 

589.68 

Missile  #41. 

Launch  time:  0  sec. 

Time  of  first  burst:  0  sec. 
Fireball  rise  not  considered. 
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TABLE  J-II 


DATA  USED  TO  PLOT  FIGURE  9 


step  ■  t 
*  max 

Time  step  ■ 

t  /2 

max' 

j 

t2(°k) 

j 

mmsam 

■n 

1 

325.38 

1 

0.25 

297.48 

2 

377.35 

2 

0.75 

325.37 

3 

446.83 

3 

1.25 

354.82 

4 

505.59 

4 

1.75 

377.30 

5 

575.00 

5 

2.25 

409.94 

6 

640.35 

6 

2.75 

446.75 

7 

697.09 

7 

3.25 

475.00 

8 

709.62 

8 

3.75 

505.45 

9 

701.93 

9 

4.25 

543.49 

10 

685.22 

10 

4.75 

574.80 

11 

660.33 

11 

5.25 

601.23 

12 

630.17 

12 

5.75 

640.16 

13 

597.60 

13 

6.25 

674.26 

14 

6.75 

696.77 

15 

7.25 

706.13 

16 

7.75 

709.25 

17 

8.25 

706.72 

18 

8.75 

701.46 

19 

9.25 

693.58 

20 

9.75 

684.81 

21 

10.25 

673.15 

22 

10.75 

659.95 

23 

11.25 

645.37 

24 

11.75 

629.86 

25 

12.25 

613.80 

26 

12.75 

597.31 

27 

13.25 

580.40 

28 

13.75 

563.43 

29 

14.25 

546.58 

30 

14.75 

529.71 

I 


Missile  #41. 

Launch  timet  0  sec. 

Time  of  first  burst:  0  sec. 
Fireball  rise  considered. 
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TABLE  J-III 


DATA  USED  TO  PLOT  FIGURE  10 
(Cumulative  Reaulta  Given  in  Table  J-II) 
(Noncumulatlve  Results) 


Burst  #1 

Burst  #2 

Burst  #3 

Burst  #4 

j 

t2(°k) 

j 

t2(°k) 

j 

T2(°K) 

j 

T2(°K) 

1 

325.38 

1 

324.33 

1 

321.99 

1 

326.51 

2 

374.31 

2 

371.75 

2 

366.19 

3 

376.54 

3 

392.41 

3 

388.05 

3 

379.86 

3 

390.88 

4 

400.50 

4 

394.30 

4 

383.78 

4 

393.82 

5 

405.42 

5 

396.80 

5 

384.74 

5 

393.34 

6 

406.96 

6 

396.31 

6 

382.91 

6 

389.89 

7 

405.14 

7 

393.12 

7 

378.67 

7 

383.91 

8 

400.62 

8 

387.40 

8 

372.42 

8 

376.05 

9 

393.78 

9 

379.77 

9 

364.52 

9 

366.51 

10 

386.06 

10 

371.68 

10 

356.54 

10 

357.00 

Miasile  #41. 

Launch  timet  0  aec. 

Time  of  firat  bur at*  0  aec. 
Fireball  rise  considered. 


Appendix  K.  Computer  Programs 

This  appendix  contains  two  computer  programs.  The 
first  one  presented  is  Program  Therm,  written  in  FORTRAN77. 
This  program  calculates  the  probability  of  survival  for 
the  cumulative  case.  Preceding  the  program  is  a  list  of 
variables  and  a  flow  chart.  The  second  program,  written 
in  Basic,  calculates  the  probability  of  survival  in  the 
noncumulative  case. 


R 

I 

I 

1 

1 

H 

i 

I 
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Variable  Lilting  for  Program  THERM 


NAME 

UNITS 

DESCRIPTION 

Main 

Program i 

a 

j/m2-°K 

>  -  epPd 

alfa 

a» 

absorptivity  of  aluminum 

alpha 

HI 

alt (11) 

- 

array  containing  altitude  at 
time  stop  j 

ang(ll) 

radians 

array  containing  missile 
flight  path  angle  at  time 
step  j 

angdata(71) 

degrees 

array  containing  data  for 
missile  flight  path  angle 

Bl,B2,B3,B4 

- 

cell  number  of  burst  1,2,3, 4 

beta 

- 

-  (l/4.108)*ln(l#k/l-g) 

Cp 

J/Kg-*K 

specific  heat  oapaoity 

d 

m 

thickness  of  missile  skin 

dCT (4 , 11) 

- 

array  containing  A  CT  for  burst 
k  at  time  step  j 

drd 

m 

missile  down-range-distance 
from  first  row  of  silos 

h(ll) 

J/m2-s-°K 

array  containing  heat  transfer 
coefficient  at  time  step  j 

Hfb (4 , 11) 

m 

array  containing  height  of 
fireball  for  burst  k  at 
midpoint  of  time  step  j 

HIT 

°K 

highest  maximum  temperature 
reached  during  calculation 
of  P 
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NAME 

UNITS 

DESCRIPTION 

** 

°K 

sure-kill  intensity 

I»8 

°K 

sure- safe  intensity 

j 

mm 

time  step 

k 

at 

burst  # 

LowT 

®K 

lowest  maximum  temperature 
reached  during  calculation 
of  P 

It 

• 

missile  launoh  time 

maxb 

- 

maximum  number  of  bursts 
(<  4) 

Mn 

- 

missile  # 

Pa 

- 

probability  of  survival 

rho 

kg/m3 

density  of  aluminum 

rhocep(lO) 

— 

array  containing  values  of  the 
distance  to  the  centroid  of 
cell  i 

aumPd 

- 

sum  of  Pd(MaxT)  for  each  burst- 
missile  configuration 

sx(100) 

m 

x-coordinate  of  silo  position 

«y (100) 

m 

y-coordinate  of  silo  position 

TeO 

°K 

ambient  air  temperature  at 
time  of  first  burst 

tbl 

a 

time  of  first  burst 

Temp (11) 

°K 

array  containing  ambient  air 
temperature  at  time  step  j 

tf 

- 

thermal  fraction 

theta (10) 

radians 

array  containing  values  of  the 
angular  Location  of  the  centroid 

cell  of  i 
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NAME 

UNITS 

DESCRIPTION 

ul 

.1 

wm 

upper  limit  to  #  of  cells 
considered  (<  I'M 

vdata(Sl) 

ft/s 

oirray  containing  da-;*  for 
missile  velocity  as  a  function 
,i  of  time 

vel(ll) 

m/s 

array  containing  data  for 
missile  velocity  at  time 
step  j 

W 

MT 

yield  in  megatons 

xdata(51) 

ft.. 

array  containing  data  for 
missile  altitude  as  a  function 
of  time 

zdata(51) 

ft 

array  containing  data,  for 
missile  distance  from' Silo 
as  a  function  Of  time 

Subroutine* j 

1 

.M  '  1 

Only  ambiguous  variables  are  listed  > 

TCalo t 

'i  •  .  ■>  ■■  ■ 

nb 

«* 

1  •  » 1 

number  of  bursts 

tb 

s 

time  of  burst 

ab 

- 

silo  of  burst 

newbst 

- 

logical  variable t  TRUE  when 
another  burst  occurs  during  j 

z 

- 

z  ■  (InMaxT  -  a')  /  &' 

C.lVar : 

dgzx,dgzy 

m 

x  and  y  coordinate  of  desig- 
nated  ground  zero 

Mny 

m 

missile  y  coordinate 
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NAME 

UNITS 

DESCRIPTION 

Inltcht 

tm 

..  ,i 

i 

a 

tim«  Into  missile  flight 

FBCalc* ,:l ,, 

tb 

a 

tima  of  burst 

tp ;■ 

- 

normalized  time 

t  ; 

a 

time  after  burat 

CTd , CTU 

' 

fraction  of  thermal  energy 
emitted  up  to  beginning  and 
up  to  end  of  time  step  j 

MsChar i 

Mnx 

m 

missile  x  coordinate 

td,t,u 

a 

integer  times  above  and  below 

t 

a 

time  of  interest 

HCalai 

op 

J/kg-#K 

specific  heat  capacity  of  air 

Re 

- 

Reynolds  number 

Pr 

•  » 

tm 

Prandtl  number 

Nu 

am 

Nusselt  number 

Ta 

°K 

ambient  air  temperature 

P 

N/m2 

ambient  air  pressure 

rhoa 

kg/m3 

ambient  air  density 

mu 

kg/m-a 

ambient  air  viscosity 

kappa 

J/m-a-°K 

ambient  air  conductivity 

xxn 

m 

point  along  missile  skin 
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Calculi  tUaaila  Charactaristica, 
h,  and  Amoiant  Air  Tamparatura 


**»#«*»«»#*>*#*!  #*#*ftftftt*»*#*««ftft*»»*»*ft#*»*rt»ftftft**»*«**ft»**ftft#»»#ft»ft»» 
ft  » 

*  PRGGRAFI  Tharm  * 

*  » 

»  Written  by  2Lt  Barbara  A.  Hall*  GNE-8BM  * 

•  Air  Fore*  Inatituti  of  Tachnology  * 

•  Fall  Quartar,  19eA  » 

•  » 

*  Thia  program  oaloulataa  tha  probability  of  survival  of  a  miaaila  * 

*  aubjactad  to  tharmal  radiation  from  a  aarias  of  bursts  for  tha  * 

*  ojmulatiua  cats.  Tha  aoanario  modaled  in  thia  protean  ia  tha  * 

*  Paacakaapar  cloaa-apaead  baaing  formation,  or  Danaa  Pack.  * 

*  Drtaila  of  tha  aoanario  ara  aa  follow# i  * 

«  ft 

*  Thraat  i  Ualk  attncls  atarting  at  ailo  fl  and  procaading  a vary  2  * 

*  atoonda  to  aueoaaaiva  ailoa  * 

ft  Surfaea  burata  * 

n  Yiald  of  aaeh  burst  i  2  fiT  * 

*  Aiming  arror  only  for  daaignatad  ground  zero  * 

«  » 

*  riiaaila  i  Alumnum  akin  * 

*  las  -  619  K,  Iak  <  809  K  * 

«  ft 

ft  Tha  program  will  prompt  tha  uaat  for  tha  following  information!  * 
a  * 

*  flissila  I  * 

*  Launch  tlma  (aao)  * 

*  Unbar  of  bursts  to  ba  oonoidarad  t  must  ba  <■  4  * 

«  Timo  of  first  burst  (aao)  i  must  ba  an  suan  numbar  and  * 

»  >«  to  tha  launch  timo  * 

«  Numbar  of  oalla  to  ba  eonaidsrad  «  batwasn  1  and  10  * 

*  Option  for  firaball  riaa  t  risa  ■  1  * 

«  ft 

*  Raaults  of  tha  calculation  will  ba  plaead  in  tha  fila  "thm.dat"  » 

*  Thm.dat  will  contain  tha  abova  information  plus  tha  lowaat  and  » 

ft  highaat  maximum  tsmpsraturaa  calculated  during  tha  run  and  tha  * 

ft  miaaila 'a  probability  of  survival.  * 

*  ft 
ftftftftftftftftftftftftftftfttftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftftft 

PROGRAfl  Tharm 

INTEGER  tb1,Mn,maxb,B1 , 02, 83, 84, ul, riaa 

REAL  I*#, Iak, Y, lii, alpha, bataftmax.tf.rhOfalfa.Cp.d.H, it 

REAL  vdatc,zdata,xdata,angdftt,sx,sy,rhocop(10) ,thata(10) 

REAL  val,alt,drd,ang,h,Tamp,TaC,dCTfHfb,LouT,HiT,sumPd,Ps 
COITION  /blockl  /vdata  ( 31 ) ,  zda  ta  ( 51 ) ,  xdata  ( 81 ) ,  angdat  ( 71 ) 

CQrtmN  /block2/ox(ina),«y(1C0) 

COITION  /block3/val(11  ),alt(11  ),drd(11  ),r,ng(11  ),h(11  ),Tamp(11 ) 
COITION  /block4/dCT(/t,11  ),Hfb(4,11 ) 

COITION  /blockd/tf ,alfa,Y,o, alpha ,bata 
DATA  81,82,03,84  /4»0/ 


I 


I 


DATA  rhocep  /0.Q,.7109,.7109, .7109, .7108,1.508, 1,509, 1.509, 
+  1.(309,1.509/ 

OATA  theta  /O.O.. 70530, 2. 3561ft, 3. 92698, 5. 60776, 0,0,1 .25633, 
4  2.51327,3.76991,5.02854/ 


C  Main  Program 
las  >  619. 

Isk  -  809. 

Y  -  2000. 

alpha  -  (1./2.)*AL0G(Ias*Iak) 
bata  -  (1./4.1O8)»AL0G(Isk/Iai) 
tmax  -  .0417*(Y*S44) 
tf  •  .18 

C  Values  for  Aluminum  i 
rho  a  2700. 
alfa  -  .50 
Cp  a  900. 

d  a  .001 

a  ■  Cv*rha»d 
CALL  Gdsta 

0PEN( 2 , fils- 1 thm.dat ' , status- 'new* , accsssa ' aaqusntial ' , forma 
+  ' formatted ' ) 

1  WRITE  (*,'(**)')  1  Mlsslls  f  (Enter  lOlor  greater  to  sxit)i  ' 
READ  (V(BN,X4)')  rtn 

IF  (l>ln  .GT.  101)  GOTO  100 

WRITF.  (•.'(A*)*)  1  Designated  launch  time  (sec)  »• 

READ  (*,'(BN,F4.1)')  It 

2  WRITE  (V(A *)')  '  Number  of  burata  fc«4)  to  be  eonsideredi  ' 
READ  (%'(&N,I1 )')  maxb 

IF  (maxb  .GT.  4)  GOTO  2 

3  WRITE  (*, * (A±) • )  '  Time  of  first  burst  (must  be  >-  It  &  evert) «  ' 

REAr;  ,'(BN,I3)')  tbl 

IF  (tbl  .LI.  It  .OR.  INT(tb1/2)  .NE.  tbl/2)  GOTO  3 

4  WRITE  (*,'(**)')  '  Number  of  cells  to  be  eonsideredi  ' 

READ  (*,'(8N,I3)')  ul 

IF  (ul  .LT.  1  .OR.  .GT.  10)  GOTO  4 

5  WRITER*,' (At)')  '  Type  1  if  went  fireball  rise,  0  if  noti  ' 
REAOtVlBN.K)')  rise 

IF  (rise  .NE.  0  .OR.  rise  .NE.  1)  GOTO  5 

CALu  InitCh (tbl  ,lt,tmax,IVi,TaO) 

W  -  Y/1000. 

CALL  FBCalc(maxb,tb1,tmax,rise,W) 

C  Lout  aiid  HIT  are  initialized 
LowT  -  I.Oe+B 
HiT  ■  0.0 
sumPd  -  0,0 

C  Calculate  ths  maximum  skin  temperature  and  probability  of  damage 
C  far  a  total  of  10**maxb  burst  -  missile  configurations. 

C  Bk  ■  current  oell  number  of  burst  k 
IF  (maxb  .EQ.  1)  THEN 
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DO  10  B1  .  1,ul 

CALL  TCalc(B1  ,B2,B3,B4,Mn,tmax,tb1  ,n*»xb, TaO,rhoceptthata, 

+  aumPd,LowT,HiT,lt) 

10  CONTINUE 

ELSE  IF  (mxb  .EQ.  2)  THEN 
DO  20  B2  •  1,ul 
DO  20  01  .  1,ul 

CALL  TCaic(81 ,82,B3,B4,Rn,tmax»tb1  ,maxb,TaD,rhocap,thatt, 
4  aumPd,LowT,HlT,lt) 

20  CONTINUE 

ELSE  IF  (maxb  .EQ.  3)  THEN 
00  30  B3  ■  1  ,ul 
00  30  B2  -  1,ul 
DO  30  B1  -  1,ul 

CALL  TCalc(S1  ,B2,B3,B4,Mn,tmax,tb1  , maxb, TaO.rhocsp, that*, 
+  iurrPd,LowT ,  HIT  .It) 

30  CONTINUE 

ELSE 

DO  40  B4  ■  1,ul 
DO  40  B3  -  1,ul 
DO  40  B2  -  1,ul 
DO  40  01  -  1,ul 

CALL  TCaic(B1  ,B2,B3,B4,Mn,tnax,tb1  ,maxb,TaQ,rhocop,thata, 
+  aunPd'LouT'HiT.lt) 

40  CONTINUE 

END  IF 


C  Writa  to  fila 
WRITE(2,89) 

88  FORMATC  ') 

URITE(2,0O)(Vi,lt,maxb,  tbl  ,ul 

90  FORMATC  «X  #\I2f'  launch  t  -  *,FB.2,'  maxb-M2,'  tbl- 

+  ,12. '  |  calls  i  ',12} 

IF  (claa  .EQ.  1)  WRITE(2,81) 

91  FORMATC  Flraball  riaa  conaidarad') 
WRITE(2,97)LowT,HlT 

97  FORMATC  Lowaat  T  -  ',F9.3,'  Highaat  T  -  '.F9.3) 

Pa  ■  1.0  -  (sumPd/(ul**maxb)) 

URITE(2,9S)Pa 
95  FORMATC  Pa  ■  '.FB.4) 

GOTO  1 


100  CLOSE  (2) 
END 


»  SUBROUTINE  TCalc  >  Temperature  Calculation  * 

u  a 

*  Ptsaad  t  81,  B2,  B3,  B4,  Mn,  tmax,  tb1t  maxb,  TaO,  rhocap(10)»  * 

«  thata  (10),  aurPd,  LowT,  HiT,  It  » 

*  COMMON  i  ax(IOO),  ay(IOO),  val(11),  alt(11),  drd(11),  ang(11),  * 

*  h(11),  Temp(11),  dCT(4,11),  Hfb(4,11),  tf,  alfa,  Y,  a,  * 

*  alpha,  bata  * 

»  a 

*  Called  by  i  Wain  program  * 

*  a 

*  Purpoaa  i  Calculataa  tha  maximum  tamparatura  and  corraaponding  * 

*  probability  of  damage  for  tha  burat  miaaila  configur-  * 

*  ation  defined  by  eella  Bk  * 

a  a 

*  Calls  i  CIVar  -  Call  Variables  * 

*  MsChar  -  Missile  characteristics  (IF  NECESSARY)  * 

a  a 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

SUBROUTINE  TCalo(B1 ,B2,B3,B4, Mn, tmax, tb1,maxb, TaO, rhocep, 

♦  theta,  aunM,  LowT,  HiT, It) 

INTEGER  B1  ,B2,B3,B4,Mn,tb1  ,nb,callno,tb,ab 

REAL  tmax, naxT,a,tf,Y, alfa, dQ,T1 ,T2,3R,CF,Tau,TaO,aunPd,Pdl,Z 

REAL  axtay,rhooap(1Q),thata(10),val,alt,drd,ang,h,Tamp,dCT,Zp,Pz 

REAL  Hfb, It, v, height, x, phi, t 

REAL  LowT, HiT 

LOGICAL  newbst 

COMMON  /blook2/ax(10Q),ay(100) 

COMMON  /block3/vel(1 1 ) ,alt(1 1 ) ,drd( 1 1 ) ,ang(1 1 ) ,h(1 1 ) ,Temp(1 1 ) 
COMMON  /block4/dCT(4,11),Hfb(4,11) 

COMMON  /blookS/tf,alfa,Y, a, alpha, beta 

nb  ■  1 
j  -  0 
T1  -  TaO 
T2  ■  T1  +  1.0 

C  Condition  for  continuing  calculation  i 
10  IF  (T2  .GE.  T1  .AND.  j  ,LE.  9)  THEN 
i  -  J  ♦  1 

C  If  new  burst  has  detonated  during  j,  increase  I  of  burste  by  1 

IF  (nb  .LT.  maxb  .AND.  J*tmax  .GE.  2*rt>)  THEN 
r<6  m  nb  +  1 
newbat  ■  .TRUE. 
a.SE 

newbat  »  .FALSE. 

END  IF 
dO  -  0.0 

C  Calculate  dQ  for  each  burat  k  at  time  step  j 

DO  20  k  -  1,nb 
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tb  -  tbl  ♦  2*(k-1 ) 
ab  -  (tb+2)/2 
IF  (k  .EQ.  1 )  Wfl 
callno  ■  81 

ELSE  IF  (k  .EQ.  2)  THEN 
callno  ■  82 

ELSE  IF  <k  .EQ.  3)  THEN 
oillno  ■  83 

ELSE  IF  (k  .EQ.  4)  THEN 
Oillno  a  B4 
END  IF 

C  If  naw  bunt  hia  Juat  occurred,  must  calculate  missile 

C  charaotariatloa  at  midpoint  between  tbk  and  J*tmax 

IF  (naubat  .AND.  k  .EQ.  nb)  THEN 

t  >  tbl-lt  +  j*tmax  -  (tb1+j*tmex-tb)/2. 

CALL  MaChar(t,w,height,xtphl,ax(Mn))  . 

CALL  ClVar(SRlCFlTaulvthaiohtlXtphitax(r>b)»ay(ab)» 

+  ay  ((In),  rhocap(  callno ) ,  theta  ( eallno )  ,Hfb(  k ,  J ) ) 

ELSE 

CALL  ClVar(SR(CFtTaulual(j),alt(j)ldrd(j),ang(j),sx(ab) 
+  ,ay(ab),ay(Cln),rhocap(oallno),thata(callno), 

+  Hfb(k,j)) 

END  IF 

dQ  -  dQ+dCT ( k , J ) *tf *Y*Tau*CF*4 .1 8Ba1 2/  ( 4.  *3 . 1 41,38*  ( SR**2 ) ) 
20  CONTINUE 

T2  -  (T1»(a-h(J)«tmax/2.)+h(j)*tmax*Tawp(d)«lfa*dQ)/ 

+  (a+h(J)»tmax/2.) 

IF  (T2  .GT.  T1)  T1  -  T2 
GOTO  10 
END  IF 
MaxT  -  T1 

IF  (MaxT  .LT.  LowT)  LowT  -  MaxT 
IF  (MaxT  .GT.  HIT)  HIT  ■  MaxT 

C  Caloulata  tha  probability  of  damaga  and  add  to  running  sum 
Z  -  (ALQG(MaxT)  -  alpha) /beta 
Zp  ■  AflS(Z) 

Pa  ■  1.-1  ./(2.*(1.+.1fl6884*Zp+.11S104*(Zp»*2)+. 000344* 

+  (Zp**3)+.01BS27*(Zp**4) )**4) 

IF  (Z  .GE.  0.0)  THEN 
Pdi  -  Pi 
ELSE 

Pdi  -  1.0  -  Pi 
END  IF 

sumPd  ■  sumPd  +  Pdi 
END 
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SUBROUTINE  CIVar  -  Call  Variablaa 

Paaaad  i  SR,  CF,  Tau,  val(j),  alt(J),  drd(J),  «ng(j),  ax(ab), 

ay(ab),  sy(Mn),  rhocap(callno),  thata(callno),  Hfb(k,J) 

Callad  by  i  Subroutina  TCale 

Purpoaa  t  Caloulataa  3R,  CF,  and  Tau  uaing  tha  givan  miaaila 
and  butat  oharaetarlatioa 


a 

a 

a 

a 

« 

a 

a 

a 

a 

a 

a 


SUBROUTINE  CIVar ( SR ,CF , Tau , u , 2 , x , phi ,dgix , dgzy ,Mny , rocapi , 

+  thataiffbh) 

REM.  SR,CF,Tau,v,z,x,phi,dgzx,dgzy,Mny,rocapi,thatai,fbh 
REAL  rhoi,xb,yb,gr,eoephi 

C  Notei  CEP  ■  200  natara 
rhoi  ■  rocapi*2Q0 
xb  ■  dgtx  +  rhoi*C03(thatxi) 
yb  ■  dgzy  4  rhoi«SIN(thatai) 
gr  a  3QR7((Mny-yb)**2  +  (x-xb)**2) 

SR  -  SQRT(gr*«2  +  (z-fbh)**2) 

IF  (v  ,EQ.  0)  THEN 
CF  a  t.O 
ELSE 

coaphi  -  ( (x-xb)*v»COS(phi)+(z-fbh)»V»SIN(pbi) )/(SR*v) 

CF  a  SQRT(1,  -  coaphi**2) 

END  IF 

TauaEXP(-.a24BB-6.A3Sa-S«3R-1 ,407a-8*(SR*»2)+1 .7B2a-14«(3R**3)) 
END 


SUBROUTINE  InItCh  -  Inlblal  Charaetariatlca 


Paasad  i  tbl,  It,  tnox,  Pin,  TaO 

COmON  i  ax(IOO),  ay (100),  vel(11),  alt(11),  drd(11),  ang(11), 
h(11),  T«p(11) 

Called  by  t  Main  program 

Purpoaa  t  Calculates  and  aboraa  mlaaila  char actariat lea,  h,  and 

ambient  air  tamparatura  for  aaoh  time  abap  up  bo  j  ■  11. 
Valuee  ara  calculabed  ab  bha  tlma  abap  midpoint  b  ■ 

(tbl -lb)  +  (j-.S)*tmax 

Calls  t  naChar  -  Mlaaila  charaebaritbica 

HCalc  -  Heat  transfer  coefficient  calculation 


SUBROUTINE  InitCh(tb1 ,lt,tmax,Mn,TaO) 

INTEGER  tbl, Fin 

REAL  lt,tmax,tO,tm,vQ,zO,xQ,phiQ,hO,TaO 
REAL  vel,alt,drd,ang,h,Tamp,ax,ay 
COWON  /biocW/ax  (1 00 ) ,  ay ( 1 00 ) 

C  Only  ax  will  bo  used  in  thie  subroutine 

COWON  /block3/vel(11),alt(11),drd(11),ang(11),h(11),Temp(11 ) 

C  Find  air  tnapsrabura  at  to 
to  •  tbl  -  It 

CALL  NaChar(bOluO,iQ,xO,phiO,ax(W)) 

CALL  HCalc(vO,zO,hO,TaO) 

C  Calculate  mlaaila  characteristics,  h,  and  Tamp  for  midpoint  of  each 
time  stop  J 

DO  10  J  -  1,11 

tm  ■  to  +  (J-.S)»tmax 

CALL  MoChsr(tm,vel(j),ait(J),drd(J),ang( J),ex(Ph)) 

CALL  HCalo(val(J),alt(J),h(J),Tamp(J)) 

10  CONTINUE 
END 
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t  • 

►  SUBROUTINE  FBCalc  -  Fireball  Calculation  * 

>  • 

1  Paaaed  i  maxb,  tbl,  tmax,  riaa,  tii  » 

-  COMMON  t  dCT(6,11),  Hfb(B,11 )  • 

»  » 

>  Called  by  i  Main  Program  * 

t  » 

1  Purpoaa  i  Caleulataa  dCT  and  Hfb  for  aach  burat  k  at  aach  tine  * 

►  atap  j  * 

t  » 

►  Calla  i  Function  CT  -  caleulataa  CT  at  tima  tp  * 

>  a 
>*#*«»*#»*»*»*a*»##a»*»#»##»a**»**#*##a»«»a#****#**»**#*#»*###»#**»**»* 


SUBROUTINE  FBCalo(maxb,tb1 , tmax, riaa, W) 

INTEGER  maxb , tbl ,nb,tb, riaa 
REAL  tmax, CT,dCT,Hfb,CTu, CTd, W,t 
LOGICAL  nawbat 

COMMON  /block4/dCT(4,11 ),Hfb(4,11 ) 

C  Inltialiia  arraya 
00  B  k  •  1,4 
OO  B  j  -  1,11 
dCT(k.j)  -  0.0 
Hfb(k,J)  -  0.0 
B  CONTINUE 
nb  ■  1 

C  For  aaeh  tima  atap  ji 
00  10  J  -  1,11 

C  If  a  new  burat  haa  datonatad  during  j,  incraaaa  I  of  burata 

IF  (nb  .LT.  maxb  .AND.  j*tmax  .QE.  2*nb)  THEN 
nb-nb+1 

nawbat  -  .TRUE. 

ELSE 

nawbat  •  .FALSE. 

END  IF 

C  For  aaoh  burat  k  that  haa  occurad  t 

00  20  k  -  1,nb 
tb  -  tbl  ♦  2*(k-1 ) 

C  Find  dCT  for  burat  k  at  tima  tp  -  t/tmax 

tp  ■  (tbl  J*tmex  -  tb)/tmax 
IF  (tp  .LE.  10)  THEN 
CTu  -  CT(tp) 
tp  ■  tp  -  1.0 
IF  (tp  »GT.  0)  THEN 
CTd  ■  CT(tp) 

ELSE 

CTd  *  0.0 
ENO  IF 
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dCT(k.J)  ■  CTu  -  CTd 
ENO  IF 

C  If  required,  find  fireball  height  for  burst  k  at  time  t 

IF  (rise  .EQ.  1)  THEN 

C  If  i  mu  burst  has  just  occur  rad  during  J,  than 

C  Hfb  is  calculated  for  tha  midpoint  of  the  time 

C  from  tbk  to  t  ■  JHmax.  Otherwise,  Hfb  is 

C  calculated  at  tha  midpoint  of  time  stop  j 

IF  (newbat  .UNO,  k  .EQ.  nb)  THEN 
t  ■  (tbl  ♦  J*tmax  -  tb)/2. 

ELSE 

t  ■  tbl  +  (J-.5):*tmax  -  tb 
ENO  IF 

Hfb(k«J)  ■  2194Q,B*(!ll**»177)*(1  »-(1  »-t/240.  )**2) 

ENO  IF 

20  CONTINUE 

10  CONTINUE 

ENO 


129 


l 


»  < 


i 


1 

I 


I 


FUNCTION  CT 


Passed  i  CT,  tp 


Called  by  i  Subroutine  FBCalo  -  Fireball  calculations 

Purpose  i  Calculates  CT  at  normalised  time  tp  from  equations 
determined  using  linear  regreaaion  on  Figure  7*64  of 
Glaeatom  and  Dolan  (Glaeatone  and  Dolan,  1977i311) 


FUNCTION  CT  (tp) 

IF  (tp  .LE.  .78)  THEN 

CT  «  -,02*tp  +  ,24*(tp**2) 

EL3E  IF  (tp  »GT«  .75  .AND.  tp  .LE.  1.5)  THEN 
CT  -  .32*tp  -  .12 

ELSE  IF  (tp  .OT.  1.5  .AND.  tp  .LE.  2.5)  THEN 
CT  ••  - . 25721 9+ . 55641  S*tp- .  0650033*  ( tp**2 ) 
ELSE  IF  (tp  .GT.  ?.5  .ANO.  tp  ,LT.  10.0)  THEN 
CT  «  .335808+, 0949904*tp-. 0049481 4# (tp«H»2) 
ELSE  IF  (tp  .Eq.  10,0)  THEN 
CT  -  o80 
END  IF 
END 


R 

1 

1 

I 

I 

I 

R 


SUBROUTINE  CtaChtr  -  Fliioilt  Characteriatica 


* 

» 

Pawed  i  t,  vel(j),  alt(j),  drd(j),  ing(J),  ax(Mn)  * 

COMMON  i  vdata(BI)*  zdata(51),  xdata(BI),  angdat(Sl)  * 

» 

Callad  by  i  Subroutine  InitCh  -  Initialize  charaoterietica  * 

# 

Pucpoau  i  Calculate*  miaaila  velocity*  altitude ,  down-range-  * 

dietance,  and  flight  path  angle  at  time  t  uaing  data  * 

fllaa  created  front  plot  of  miaaila  behavior  (Figure  1 )  * 


SUBROUTINE  MeChar(t,v,*,x,phi,Mnx) 


REAL  t,v,z,x,phi,mnx 
INTEGER  td,tu 

HEM.  udata . sdata . xdata . anode t 

COMMON  /blcekl /vdata( 81 ) , sdata ( 51 ) , xdata( 51 ) ,angdat(71 ) 

C  Determine  lower  and  upper  limit  of  interpolation 
td  -  INT(t) 
tu  ■  td  ¥  1 

C  Find  velocity*  altitude*  and  diatanea  from  ailo 
If  (t  .GT.  SO)  THEN 
v  >  108»*t  -  1350* 
x  m  24B0*«t  -  55000. 
x  «  3700. *t  -  127000. 

ELSE 

C  Uaa  linear  interpolation  to  calculate  value 

v  ■  ( t-td ) * ( udata ( tu+1 ) -udata ( td+1 ) )  4  vdata(tdfl) 
z  ■  (t~td)»(zdata(tU4l)-zdata(td4-1))  4-  zdata(td4-1) 
x  ■  (t~td)M(xdata(tu4l)>xdata(td4-1))  4  xdata(td4l) 

END  IF 

C  Find  flight  path  angle 
IF  (t  .GT.  70)  THEN 
phi  >  30. B  -  ,1*(t-70) 

ELSE 

phi»(t-td)*(angdat(tu+1  )-angdat(td-t-1  ))+«ngdat(td4l ) 

END  IF 

C  Convrrt  to  metera  and  radlane  and  calculata  down-range-diatanca  from  x«0 
v  -  u*.304S 
z  «  z*.3O40 
x  ■  x*.3O40  4-  Mnx 
phi  -  phi*3.14138/1 80. 

END 


SUBROUTINE  HCalc  -  Heat  Tranafac  Coefficient  Caloulation 

Paaaad  i  vel(J),  alt(j),  h(J),  Tamp(J) 

Callad  by  i  Subroutirw  InitCh  -  Initializa  charaotsriatioa 

Purpose  i  Calculates  h  using  missile  velocity  and  ambient  air 
properties  at  miaaila  altitude.  Air  propartiaa  ara 
found  using  aquations  from  US  Standard  Atmosphere 
(NOAA,  1878 i 0-30) 


SUBROUTINE  HCalc(v,2,htc,Ti) 

REAL  v.z.hto.Ts 

REAL  opJftu»Pr*Nu,xm*P»rhoaimu»kappa 

C  Find  ambiant  air  conditions  at  z  using  US  Standard  Atmoaphara  aquations 
IF  (i  .LT.  11000.)  THEN 
Ta  ■  2S8.1S  -  .006848*2 
P  -  101 300. *((2B8.15/Ta)**(-. 0341 8A/.OOQS45)) 

ELSE  IF  (2  .GE.  11000.  .AND.  2  .LT.  20000.)  THEN 
Ta  •  216.68 

P  -  22880.*EXP(-.034184»(2-1100Q.)/218.B8) 

ELSE  IF  (z  .GE.  20000.  .AND.  2  .LT.  32000.)  THEN 
Ta  *  218,88  ♦  .001*(z-20000.) 

P  *  5628, *( (216. SS/Ta)**( ,034164/«(J01 ) ) 

ELSE  IF  (z  .GE.  3SOOO.  .AND.  z  .LT.  47000.)  THEN 
Ta  «  226.66  4  ,Q02b»(z-3200C.) 

P  >  088.8»((22a.88/Ta)»*(. 0341 64/ .0028)) 

ELSE 

WRITE(*,*)**  >  47000' 

END  IF 

rhoa  -  .003404*P/Ta 

mu  •  1.45Ba-8*(Ta««1.3)/(Ta  ♦  110,4) 

kappa  ■  2.84B3Be-3*(Ta**1 ,8)/(Ta  ♦  24S.4*(10**(-12./Ta))) 

C  Calculata  tha  hast  tranafar  ooafficiant  at  point  xm 
op  -  240, *4.1 04 
xm  x  S.S 

Ra  x  rhoa*v*xm/iiiu 

Pr  r.  mu*cp/kappa 

IF  (Ra  .LE.  500000.)  THEN 

Nu  ■  .332*(Pr**(1./3.))*(lia**.5) 

ELSE 

Nu  •  ,0296»(Pr**(1./J.))*(Re**(4./5.)) 

ENO  IF 

htc;  *  Nu*kappa/Xffl 
ENO 
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SUBROUTINE  GOate  -  Gat  Data 


COWON  t  vdata(SI),  zdata(51),  xdeta(SI),  angdat(BI),  ex(IOO), 

■y(ioo) 

Called  by  i  Main  program 

Purpoaa  i  Inltlallzaa  mlaailr  oharaotarlatie  arraya  by  reading 
data  from  external  filaa  and  initiellxea  alio  position 

arraya 


SUBROUTINE  Gdata 

REN.  vdata,zdata,xdata,angdat|axlay,dx,dyirhooaptthata 
CQWON  /blockl  /vdata  ( 51 ) ,  zdata  ( 51 ) ,  xdata  ( 51 ) ,  angdat  ( 71 ) 
COWON  /block2/ax( 1 00 ) , ay ( 1 00 ) 

C  Read  In  data  from  filaa 
10  FOWAT(8N,F0.1) 

20  F0RF!AT(BN,F8.2) 

OPEN  (B.filea’veldata.txtSatatueaW, 

+  aooaaaB'aaquantial' tfarm*' formatted1 ) 

REWIND  8 
DO  30  1  -  1,81 

READ(8,1C)vdate(i) 

30  CONTINUE 

CUWg  (g) 

OPEN  ( 8 , f ila» ' altdata . txt ' .atatus-'old' , 

♦  aocaaa«'aaquantial' , form* 'format tad') 

REWIND  9 
DO  40  1  »  1 ,81 

READ(B,10)zdata(i) 

40  CONTINUE 
CLOSE  (8) 

OPEN  ( 1 0 , f ile« ’ drddata . txt ' , ita tua- ' old ’ , 

+  aooaaa» ' aaquant ial ' , forma ' formatted ' ) 

REWIND  10 
00  SO  i  -  1 ,81 

REA0(10,1Q)xdeta(i) 

SO  CONTINUE 
CLOSE  (10) 

(PEN  ( 1 1 , filaa ' dagdata • txt ' , atatnaa ' old ' , 

+  aceaeax ' aaqmmtial ' , f arm* ' formatted ' ) 

REWINO  11 
DO  BO  1  -  1,71 

REA0(11,20)angdat(i) 

SO  CONTINUE 
CLOSE  (11) 
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C  Calculate  and  atora  alio  x  and  y  poaltion 
dx  -  519.82 
dy  ■  300. 

00  70  1  -  0,85,5 
ay(l+1)  ■  0*0 
ay(i+2)  »  2.#dy 
ay (1+3)  ■  4,»dy 
•y(l+4)  •  dy 
ay(l+5)  ■  3,*dy 
70  CONTINUE 
J  -  0 

DO  80  1  •  0,38,2 
ax(J+1)  •  i*dx 
ax(J+2)  •  l*dx 
ax(J+3)  -  i*dx 
ax( J+4)  -  (1+1 )»dx 
ax(J+5)  -  (1+1 )*dx 
J  -  J  +  8 
00  CONTINUE 
ENO 
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10  REM  »***»**•  THERMAL  INTERACTION  /  Non-Cumulative  Effect  *»»**»*** 

20  DIM  WDATA(51),ZDATA(S1),XDATA{51),ANOOATA(?1),«.T(100)»SX(10Q)t3Y(100) 

30  G0SU6  1100  i  REM  Initialize  vel,alt,drd,ang,ndt,sxfand  ay  arrays 
40  PI  -  3*14159  i  DIM  CT(11) 

SO  DATA  O.Q, 0.2, 0*488, 0.5783, 0.5386,0.8871, 0.7277, 0,7884,0.7782, 0.7902, 0.80 
60  FOR  I  -  0  TO  10  i  READ  CT(I)  l  NEXT  I 
70  DIM  RHOCEP(IO) . THETA(IO) 

80  GOSUB  1020  i  REM  Call  subroutine  to  inltiallzs  N  osll  cantroid  co-orda 

90  INPUT"Yiald  (KT)"|Y  t  INPUT  %maxb,riBa"|N,MAXB,RI3E 

100  H08  -  0  t  153  •  619  I  ISK  -  809  i  CEP-  200 

110  ALPHA  -  (1 1/21  )aLQG(I93*!SK)s  SETA  -  (1  l/4.10B)«H.0G(I3K/I3S) 

120  7MAX  ■  ,0417*Y\44 

130  RH0-2700  i  ALFA-,5  I  CV-900  t  0  -.001  I REM  values  fur  A1 

140  A  ■  CV*RHQ*D  l  C  ■  ALFA«PMAX/{A*4!»PI)  l  TF  -  .18  I  XM  -  3.5 

ISO  INPUTH3ilo  I  and  daaignatad  launch  time"|SN,LT 

180  INPUTnTiM  of  first  burst  (must  be  >■  launch  time  and  avan)n|TB1 

170  LPRINTiLPRINT"Missile  #"|5N|"  Launch  t»"|LT»"  First  burst  at  tins  »"|TB1 

180  LPRINT"Y(KT)«"|Y,"Ni"|N,"riiaxb»w|flAXB,"riaa»"tRISE 

188  LOUT  •  10000001  »  HIT  •  0! 

190  GOSUB  230  t  REM  Find  cumulative  Pa 
200  GOTO  90 


220  REPI  WHHHHMHHHUHHHHHHHHMHMUHHHHHHUHHUHHHHHHUMUHHWWHMUHHHUHHHHHHUMHHMUHt 

230  RC1  Subroutine  to  calculate  cum  Ha  for  first  four  bursts  affecting  ON 

24J  FOR  NB  -  1  TO  PtAXB 

250  TB  ■  T81+2*(NB-1 )  t  SB  ■  (TB+2)/2 

280  DGZX  «  SX(SR)  I  OCZY  -  SY(SB)  i  SUPPO  -  01 

270  TO-TB  -  LT 

280  FOR  I  -  1  TO  N 

290  PRINT  MB, I 

300  3  -  0  i  T  »  TO  »  COSUB  510  I  Z  *  ALT  l  GOSUB  880  j  T1  ■  TA 

310  RtJ»1 

320  3-3  +  1  iT  - TO  +  (3-.S)*Tf"AX 

325  IF  RISE  ■  1  THEN  HFB  ■  21B40.B«(Y/1Q00l)M77«(11-(1  l-T/2401)*2) 

330  GOSUB  510  t  GOSUB  580  s  GOSUB  790  I  TEMP  ■  TA 

335  OCT  ■  CT(3)-CT(J-1) 

340  DQ  »  0CT»TF»Y»TAU*CF»4.1BBE+12/(4*PI*3R,2) 

350  T2  ■  (T1*(A  -H*TPlAX/2 )  +H*TMAX« TEMP+ALFA*DQ)/(A+H*TflAX/2) 

360  LPRINT  USING "|#  "»3M  LPRINT  USING**#####.##  "|T2 

370  REW  IF  T1  >  T2  THEN  GOTO  395 

380  T1  ■  T2  i  IF  3  <  10  THEN  GOTO  320 

390  REM 

395  IF  T1  <  LOUT  THEN  LOUT  -  T1  I  IF  T1  >  HIT  THEN  HIT  -  T1 
400  B  -  (L0G(T1 )-ALPHA)/0ETA 
410  BP  -  ABS(8) 

420  PB  -  1-1 l/(2l*(1 l+.1Q6854*BP+.115194*flP,2+,000344*BP,3+.019527*BP,4),4) 

430  IF  B>-  0  THEN  POI  ■  PB  ELSE  POI  ■  1-PB 
440  5UPP0  »  SUPPO  +  POI 

450  NEXT  I 

460  P0-(1  l/N)»3i*P0  i  PS-1 1-PO  i  IF  NB-1  THEN  CUPPS-P3  ELSE  CUPPS  -  OJPPS*PS 
460  NEXT  NB 

482  REPI  LPRINTlow  T  ■  *|iL 

485  REPI  LPRINT"Cumulatlvs  Ps  for  "iiLPRINT  USING "ff  "|FIAXB|iLPR3CNT"burets  1st 

"iClIPS 

490  RETURN 


136 


xvuux: 


ajsuiwxxsa 


§!§gg  gggsggsgsIS  §§SS§l§881§iiB 


600  HEW  *******##**«#**»*»#»»*#**»»*«H»»»«*«»**«»*»»»»*»**»»**»*»*»»»*«»**«* 
510  REM  Subrout ina  to  find  misaila  ual.alt.drd.  and  ang  at  tima  T 
520  TO  -  FIX(T)  «  TU  -  TO  +  1 
IF  T  >  50  THEN  GOTO  810 

IF  T«TD  THEN  VL-VOATA(T)iAL-ZUATfl(T)«DR-XOATR(T)«GQTO  820 
ROT 

VL  ■  ((T-TO)/(TU-TO))»(VOATA(TU)-V0ATA(TD))  +  VDATA(TD) 

AL  -  ( (T-T0)/(TU-TD) )*(ZDATA(TU)-ZDATA(TO) )  +  ZOATA(TO) 

OR  -  ((T-T0)/(TU-T0))«(XDATA(TU)-X0ATA(T0))  +  XDATA(TD) 

GOTO  820 
ROT 

VL  «  1 05*T-1 350  i  AL  «  24*0*T  -  850001  I  DR-3700*T-1 27000 t 
IF  T>70  THEN  ANG»  30.8  -  .1*(T-70)  l  GOTO  850 
IF  T  ■  TO  THEN  ANG  -  ANGDATA(T)  I  GOTO  850 
ANG-{ ( (T-TO)/(TU-TD) )»( ANGOATA(TU)-ANGDATA(TO) )+ANBDATA(TO>) 

VEL  -  VL«.3048iALT  -  AL*.3O40iDRD  -  0R*,3O4B+SX(3N)«  ANU  -  AN0*2 I *PI/360! 
RETURN 


RQf|  *MMW«#*»««*a#*####«**M*##*»»***»#**********»*#»*»»*»«******«»»*«»*  » 
ROT  Subroutlna  to  oaloulata  SR#  CFr  and  TAU  givan  Al.T 
DELTAZ  -  ALT-HFB  t  RHOI  -  RHOCEP(l)*CEP 

XB  ■  OGZX  +  RHQI*COS(THETA(l))  l  Y8  ■  DGZY  +  RHOI*SIN(THETA(I)) 

GR  -  SGR((SY(3N)-YB)*2  4  (DRD-XB)*2) 

SR  -  SQR(GR*2  +  DELTAZ'2)  I  SRX  -  DRO-XB  I  SRZ  -  DELTAZ 

IF  VEL  *  0  THEN  CF  «  1  *  GOTO  780 

COSPHI  -  (3RX*VEL*C0S( ANG)+3RZ*VELW3IN(ANG) )/(SR*VEL) 

CF  -  SQR(1 I-COSPHI'2) 

T  AU«EXP ( - .02455-6 . 43BE-05*SR-1 .407E-09*SR*2+1 ,792E-14*SR*3) 

RETURN 


REfl)  **IMUMI*#«****»*<MM**#*#«  #«*«#*#*»##«*##»*#******»*»»**«»»»*»******** 

ROT  Subroutlna  to  calculate  haat  tranafar  eoafficiant,  h  (J/ra2-*-K) 
z  ■  ALT  i  G03UB  880  I  ROT  Find  ambiant  air  rhoa,TA,mu,ka  at  z 
CP  «  240*4.184 
RE  -  RHOA*VEL*XW/nU 
830  PR  -  l*U*CP/KA 

840  IF  RE<-5000001  THEN  NU-.332*PR*(1/3)«RE,,5  ELSE  NU-.0288*PR,(1/3)*RE,(4/5) 
850  H-  NU*KA/XM 
880  RETURN 
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870  REM  »#t*»»**<H**«*««*##»M#»M*##»***##MM#«***»#*»#*##*(»#»*t******##»*» 
880  REM  Subroutine  »  US  Standard  Atmoapher*  to  47  km 
MO  IK  Z  <  11000  THEN  LK  -  -.00BS4S  i  PK  ■  101300!  I  TK  -  288.15  t  ZK  >  0 
900  IF  Z>«1 1 000  AND  Z<20000  THEN  LK«0  t  PK -22690 I  i  TK-21B.88  l  ZK-11000 
910  IF  Z>-20000  AND  Z<32Q00  THEN  LK-,001  I  PK-5528  ,  TK-216.85  i  ZK-20000 

920  IF  Z>-32000  ANO  Z<47000!  THEN  LK-.002B  iPK-888.8  t  TK-228.65  iZK-32000 

930  IF  Z  >-  470001  THEN  PRINT  "a  >-  47000,  ae  consult  NCAA  for  values" 

940  IF  LK  •  0  THEN  GOTO  950  ELSE  GOTO  980 

950  P  «  PK«EXP(-.Q341B4*(Z  -  ZK)/TK)  i  TA  -  TK  »  GOTO  970 

960  TA  ■  TK  +  LK«(Z  -  ZK)  I  P  -  PK*(TK/TA)*(.Q34184/LK)  I  GOTO  970 

970  RHOA  -  .003484«P/TA 

060  mu  >  1.45eE-06»TAa1.5/(TA  +  110.4)  t  REN  (kg/n-a) 

990  KA  -  2.84838E-03»TA*(1 ,5)/(TA  ♦  245,4*10,(-12!/TA))  l  REN  3/(m-e-K) 

1000  RETURN 


1010  REN  ***#*•»*«*»»****«****#«*+***»****»**«»*»*•»*«»********«**«**»»*#** 
1020  REN  Subroutine  to  initialize  coordinate*  of  N  equal  prob  call  grid 
1030  RH0CEP(1)-0t  I  FOR  1-2  TO  5  *  dH0CfP(I)-.7109  >  NEXT  I 
1040  FOR  I  ■  S  TO  10  l  RHOCEP(l)  -  1.809  l  NEXT  £ 

1050  THETA(1)-0I  »  THETA(2)-(Wl/4)/2  I  THETA(B)  -  01 
1080  FOR  I  -  3  10  5  l  THETA(I)  -  THETA(I-1 )+(2*PI/4)  I  NEXT  I 
1070  FOR  I  -  7  TM a  «  THETA(l)  »  THETA(I-1)+(2«PI/S)  iNEXT  I 
1080  RETURN 


1100  REN  Subroutine  to  initialize  veltalt,drd,ang,mlt,exfand  ay  array* 

1110  OPEN  "I", ft,' "Biveldata.txt" 

1120  FOR  I  -  0  TO  50  •  INPUT  f1,V0ATA(X)  I  NEXT  I  I  CLOSE  |1 
1130  OPEN  "I",|2,"Bialtdata.txt" 

1140  FOR  I  -  0  TO  50  i  INPUT  #2,ZDATA(I)  t  fOT  I  I  CLOSE  |2 
1150  OPEN  "l\|3, "Bidrddata.txt" 

1180  FOR  I  -  0  TO  50  I  INPUT  |3,XDATA(I)  t  NEXT  I  t  CLOSE  |3 
1170  OPEN  "I",|1, "Bidegdata.txt" 

1180  FOR  I  >  0  TO  70  i  INPUT  fl ,ANGDATA(I)  I  NEXT  I  »  CLOSE  |1 
1190  OPEN  "I",|1 , "BiNLO.txt" 

1200  FOR  I  -  1  TO  100  i  INPUT  |1,NLT(I)  t  NEXT  I  I  CLOSE  fl 
1210  OX  •  519.82  i  OY  -  300 
1220  FOR  I  -  0  TO  95  STEP  5 

1230  3Y(I+1 )-0  i  SY(I+2)-2»DY  i  SY(I+3)-4*0Y  t  SY(I+4)-0Y  l  SY(I+5)-3»QY 

1240  NEXT  I 
1250  3  -  0 

1280  FOR  I  -  0  TO  38  STEP  2 

1270  3X(3+1 )-I»DX«SX(3+2)-I*DXlSX(3+3)-I»DXlSX(3+4)-(I+1 )*0Xl3X(3+5)-(I+1  )»DX 
1280  3  -  3+5 
1290  NEXT  I 
1300  RETURN 
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Appendix  L.  Results  of  Dust  Shie ldlnq  Analysis 

The  following  table  presents  data  which  shows 
that  the  geometric  treatment  of  dust  shielding  does  not 
affect  the  transmittance  to  any  missiles  of  interest 
(i.e.,  missiles  past  #31).  The  values  of  <j>  given  are  the 
maximum  values  reached  during  a  four-burst,  1-cell  CEP 
scenario. 


TABLE  L-I 

RESULTS  OF  DUST  SHIELDING  ANALYSIS 


Missile  launch  timei  0 
Time  of  first  burst i  0 
1-cell  CEP 

Fireball  rise  considered 

sec 

sac 

>  20*s  T 

m 

0 

Maximum  s>  (°) 

for  Burst  # 

Missile  # 

1 

.2 

4 

11 

27.4 

18.9 

9.5 

1.5 

14 

22.9 

16.1 

8.4 

1.2 

15 

22.1 

16.1 

8.7 

1.1 

16 

19.9 

13.6 

7.0 

1.0 

21 

15.5 

10,6 

5.5 

0.7 

31 

10.7 

7.3 

3.9 

0.5 

41 

8.2 

5.6 

2.9 

0.3 

51 

6.8 

4.5 

2.4 

0.3 
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